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ABSTRACT 

Context. The Kilo-Degree Survey (KiDS) is an optical wide-held imaging survey carried out with the VLT Survey Telescope and 
the OmegaCAM camera. KiDS will image 1500 square degrees in four Liters (ugri), and together with its near-infrared counterpart 
VIKING will produce deep photometry in nine bands. Designed for weak lensing shape and photometric redshift measurements, its 
core science driver is mapping the large-scale matter distribution in the Universe back to a redshift of ~0.5. Secondary science cases 
include galaxy evolution. Milky Way structure, and the detection of high-redshift clusters and quasars. 

Aims. KiDS is an ESO Public Survey and dedicated to serving the astronomical community with high-quality data products derived 
from the survey data. Public data releases, the hrst two of which are presented here, are crucial for enabling independent conhrmation 
of the survey’s scientihc value. The achieved data quality and initial scientihc utilization are reviewed in order to validate the survey 
data. 

Methods. A dedicated pipeline and data management system based on Astro-WISE, combined with newly developed masking and 
source classihcation tools, is used for the production of the data products described here. Science projects based on these data products 
and preliminary results are outlined. 

Results. For 148 survey tiles (^160 sq.deg.) stacked ugri images have been released, accompanied by weight maps, masks, source 
lists, and a multi-band source catalogue. Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5cr in a 2"aperture) in ugri, respec¬ 
tively, and the typical r-band PSF size is less than 0.7". The photometry prior to global homogenization is stable at the ~2% (4%) 
level in gri (u) with some outliers due to non-photometric conditions, while the astrometry shows a typical 2-D RMS of 0.03". Early 
scientihc results include the detection of nine high-z QSOs, hfteen candidate strong gravitational lenses, high-quality photometric 
redshifts and structural parameters for hundreds of thousands of galaxies. 

Key words, methods: observational - astronomical data bases: surveys - galaxies: general - cosmology: large-scale structure of 
Universe 


1. Introduction 

Sensitive, wide-held astronomical surveys have proven to be ex¬ 
tremely useful scientihc resources, not only for the specihc sci¬ 
entihc use cases they are designed for, but also due to their legacy 
value and as source-hnders for the largest telescopes. With the 
arrival of VISTA in 2010 and the VLT Survey Telescope (VST) 


in 2011, the European astronomical community now has access 
to dedicated survey telescopes both in the infrared and the opti¬ 
cal. During the hrst years of operations the majority of observ¬ 
ing time on both telescopes is dedicated to a number of pub¬ 
lic surveys, large observational programs selected by ESO that 
serve the astronomical community with regular public data re- 
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lease^ On VST these are the VST Pho tometric H-a Sur vey of 
the Southern Galactic Plane (VPHAS+, [Drew et al.|2014| ), VST 
ATLAS ( [Shanks et al.|2013| ), which covers 5000 square degrees 
in the Southern Galactic Cap to similar d epth as the Sloan Dig¬ 
ital Sky Survey (SPSS, [Ahn et al.||2012| ), and the Kilo-Degree 
Survey (KiDS, jde Jong et al.||2013| ). The Kilo-Degree Survey 
(KiDS) is the largest public survey on the VST in terms of ob¬ 
serving time and will image 1500 square degrees of extragalactic 
sky in four filters, u, g, r, and i. Together with the VISTA Kilo- 
Degree Infrared Galaxy Survey (VIKING, [Edge et al.|2013| ), a 
sister survey on VISTA that observes the same area in Z, F, /, 
H and Ks , this will result in a deep nine-band data set with, in 
total, 185 minutes of observing time per survey tile over all fil¬ 
ters. Part of a long heritage of ever improving wide-field optical 
sky surveys, the combination of superb image quality with wide 
wavelength coverage provided by KiDS/VIKING will be unique 
for surveys of this size and depth for the foreseeable future. 

The main science driver for KiDS is mapping of the large- 
scale mass distribution in the Universe and constraining its ex¬ 
pansion history by means of weak gravitational lensing and pho¬ 
tometric redshift measurements. These goals put stringent re¬ 
quirements on image quality and stability, photometric depth 
and calibration accuracy. However, this high data quality will be 
beneficial for many other science cases, and secondary science 
drivers for KiDS include a varied set of topics, such as the evo¬ 
lution of clusters and galaxies, galaxy scaling relations. Galactic 
stellar halo structure, and searches for rare objects such as strong 
gravitational lenses and high-redshift quasars. How KiDS will 
benefit the pri mary and secondary s cience drivers is discussed in 
more detail in jde Jong et al.| ( [20T3] ). 

KiDS data will be released to the astronomical community in 
yearly Public survey releases of those survey tiles that have been 
observed in all four filters. In this paper we describe the data 
products of the first two such KiDS data releases (DR), covering 
a total of 148 survey tiles. Apart from the released data prod¬ 
ucts (Sect. 13, observational set-up (Sect.f^, the data processing 
pipeline (Sect. and data quality (Sect. we also review the 
scientific research that is currently being carried out by the KiDS 
team using the data presented here in Sect. Finally, the paper 
is closed with a brief summary and outlook in Sect. [ 7 ] 


2. Data description 

2.1. General properties 

KiDS is designed for weak lensing shape and photometric red- 
shift measurements, with weak gravitational lensing tomography 
as its ultimate science goal. This goal is efficiently reached at the 
VST with deep, good-seeing r-band observations taken under the 
best atmospheric conditions, supplemented with somewhat shal¬ 
lower and less sharp u, g, i images. Galaxy shapes for the lensing 
studies can then be measured on the r band data, while the other 
bands yield colours that can be used to derive photometric red- 
shifts for these same galaxies. The detailed implementation of 
these requirements is described below in Sect. typical values 
for the resolution (PSF FWHM) and limiting magnitude (5-(T 
AB magnitude in a 2"aperture) are (1.0", 0.8", 0.65", 0.85") 
and (24.3, 25.1, 24.9, 23.8) in w, g, r and i bands, respectively 
(see Sect. for a detailed description). The i band data are sig¬ 
nificantly less uniform in quality than the other bands because 
they are mostly taken when the moon is above the horizon. 


^ http: //WWW .eso.org/sci/observing/PublicSurveys.html 
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2.2. OmegaCAM and VST 

Designed to provide superb and uniform image quality over 
the full field-of-view (FOV), the combination of VST and the 
OmegaCAM wide-field imager is ideal for a survey such as 
KiDS. The VST, an alt-az mounted modified Ritchey-Cretien 
telescope, uses a two-lens field corrector and has active primary 
and secondary mirrors. OmegaCAM is a 268 Megapixel CCD 
imaging camera with a one square degree FOV. Located at the 
Cassegrain focus, it is the only instrument on the VST. The fil¬ 
ters and dewar window are part of the total optical design, with 
the latter acting as a field lens. The science array consists of 
32 thinned, low-noise 2kx4k E2V devices. With little aberra¬ 
tion over the full FOV the I5p pixel size translates to a constant 
pixel scale of 0.21" with small gaps of 25"and 85"between the 
CCDs. Wavefront sensing, using two auxiliary CCDs on either 
side of the science array, is employed to monitor and optimize 
the optical set-up in real time. Auto-guiding of tracking and ro¬ 
tation makes use of two more auxiliary CCDs. Combined, the 
VST and OmegaCAM provide a PSF equal to the atmospheric 
seeing over the full FOV down to 0 .6". More details on the VST 
and OmegaCAM can be found in jCapaccioli et al.j p012| ) and 
|Kuijken| ( |2Qll| ), respectively. 

2.3. Sky distribution 

KiDS data releases consist of ~1 square degree tiles that have 
been successfully observed in all four survey filters (u,g,r,i). To 
maximize the synergy with existing spectroscopic surveys and 
between KiDS and VIKING early on, the latter surveys both 
prioritize the observations in the Galaxy And Mass Assembly 
(GAMA, [Driver et al.|2QTT] ) fields G9, G12, G15 and G23. 

The first public data release of KiDS (KiDS-ESO-DRl) was 
issued in July 2013 and contains imaging data, masks and single¬ 
band source lists for all tiles observed in all four filters during the 
first year of regular operations (15 October 2011 to 31 Septem¬ 
ber 2012), including data taken during Early Science Time (EST, 
13 August to 15 October 2011), a total of 50 tiles. The second 
data release (KiDS-ESO-DR2) was available in February 2015 
and contains the same data products for all tiles for which ob¬ 
servations were completed during the second year of regular op¬ 
erations (1 October 2012 to 31 September 2013), a total of 98 
tiles. Since the processing pipeline used for KiDS-ESO-DR2 is 
practically identical to that used for KiDS-ESO-DRl, the tiles 
released in the former complement the data set of the latter, 
making KiDS-ESO-DR2 an incremental release. Apart from the 
data products mentioned above, KiDS-ESO-DR2 also provides 
a multi-band source catalogue based on the combined set of 148 
tiles released in the first two data releases. 

Many other fields have been observed in a subset of the fil¬ 
ters and will be included in future releases once their wavelength 
coverage is complete. Figureshows the sky distribution of the 
tiles included in KiDS-ESO-DRl/2 within the KiDS fields. A 
complete list of all tiles with data quality parameters can be 
found on the KiDS websit^El 

2.4. Data products 

For every tile in the first two data releases of KiDS the following 
data products are included for each of the bands, u, g, r, and i: 

- astrometrically and photometrically calibrated, stacked im¬ 
ages (“coadds”) 

^ http://kids.strw.leidenuniv.nl/DR2 
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KiDS-North 



Fig. 1. Sky distribution of survey tiles released in KiDS-ESO-DRl (blue) and KiDS-ESO-DR2 (green), with the full KiDS area is shown in 
grey. Top: KiDS-North. Bottom: KiDS- South. The multi-band source catalogue covers the combined area (blue -h green). The black dashed lines 
delineate the locations of the GAMA fields. 


- weight maps 

- flag maps (“masks”) that flag saturated pixels, reflection ha¬ 
los, read-out spikes, etc. 

- single-band source lists 

The multi-band source catalogue encompassing the com¬ 
bined area of the two data releases, although forming one large 
catalogue, is stored in multiple flies, namely one per survey tile. 


2.4.1. Coadded image units and gain 

The flnal calibrated, coadded images have a uniform pixel scale 
of 0.2 arcsec. The pixel units are fluxes relative to the flux cor¬ 
responding to magnitude = 0. This means that the magnitude m 
corresponding to a pixel value / is given by: 


-2.51ogio/. 


( 1 ) 


The gain varies slightly over the fleld-of-view because of the 


photometric homogenization procedure described in Sect. 

An average effective gain is provided onlin^ An exarmle of a 
FITS header of a coadded image is also provided onlim " 


2.4.2. Single-band source list contents 

For each tile single-band source lists are provided for each of 
the survey Alters. To increase the usefulness and versatility of 
these source lists, an extensive set of magnitude and shape pa¬ 
rameters are included, including a large number (27) of aperture 
magnitudes. The latter allows users to use interpolation methods 
(e.g. “curve of growth”) to derive their own aperture corrections 
or total magnitudes. Also provided is a star-galaxy separation 

^ http://kids.strw.leidenuniv.nl/DR2/data_table.php 
http: //kids. strw. leidenuniv. nl/DR2/example_ 
imageheader.txt 


parameter and information on the mask regions (see Sect. 4.4) 
that might affect individual source measurements. Details on the 
production of these source lists, including source detection and 
other measurements are discussed in Sect. 14.5.11 

Table ] A. 1 [ lists the columns that are present in the single-band 
source lists provided in KiDS-ESO-DRl/2. Note that of the 27 
aperture flux columns only the ones for the smallest aperture (2 
pixels, or 0.4" diameter) and the largest aperture (200 pixels, or 
40" diameter) are listed. 


2.4.3. Multi-band catalogue contents 


KiDS-ESO-DR2 includes a multi-band source catalogue. This 
catalogue is based on source detection in the r-band images. 
While magnitudes are measured for all of the Alters, the star- 
galaxy separation, positional, and shape parameters are based on 
the r-band data. The choice of r-band is motivated by the fact 
that it typically has the best image quality and thus provides the 
most reliable source positions and shapes. Seeing differences be¬ 
tween observations in the different Alters are mitigated by the 
inclusion of aperture corrections. Details on the production of 
this catalogue are discussed Sect. |4.5.2| and Table jB.lj lists the 
columns present in the multi-band source lists in KiDS-ESO- 
DR2. 


2.5. Colour terms 

The photometric calibration provided in KiDS-ESO-DRl and 
KiDS-ESO-DR2 is in AB magnitudes in the instrumental sys¬ 
tem. Colour-terms have been calculated with respect to the SDSS 
photometric system. 

Aperture-corrected magnitudes tak en from the multi- band 
catalogue were matched to SDSS DR8 ( jAihara et al.|2011| ) PSF 
magnitudes of point-like sources. For each filter, the median off¬ 
set to SDSS is first subtracted, rejecting tiles where this offset 
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Fig. 2. Derivation of colour-terms of KiDS photometry with respect to 
the SDSS photometric system. The distribution of stars is shown as a 
2-D histogram, overlayed with isodensity contours (green) and fit by a 
linear relation (red line). For u the contours correspond to 50, 100, 200 
and 350 stars per bin, and for g, r and i to 100, 350, 600 and 1000 stars 
per bin. The absolute calibration between KiDS and SDSS (horizontal 
offset) is arbitrary. 


exceeds 0.1 mag in any of the bands (10 tiles), in order to pre¬ 
vent poor photometric calibration to affect the results. The fit is 
performed on all points from the remaining tiles. Figureshows 
the photometric comparison and the following equations give the 
resulting colour terms: 


^KiDS - ^SDSS = (-0.050 + 0.002)(t/sDSS - ^SDSs), (2) 

^KiDS - ^SDSS = (-0.052 + 0.002)(gsDSS - ^SDSs), (3) 

^KiDS - ^SDSS = (-0.033 + 0.002)(gsDSS - ^SDSs), (4) 

iKiDS - isDss = (+0.012 + 0.002)(rsDss “ ^’sdss)- (5) 


2.6. Data access 

Data from the first two KiDS data releases can be accessed in a 
number of different ways: through the ESO science archive, from 
the Astro- WISE archive, or via a web-based synoptic table on 

the KiDS websit43 

2.6.1. ESO archive 

As an ESO public survey, the KiDS data releases are distributed 
via the ESO Science Archive Eacilit}]^ Using the Phase 3 query 
forms users can find and download data products such as the 
stacked images and source lists. 

The naming convention used for all data product files in the 
ESO archive is the following: 

KiDS_DRV.V_R.R_D.D_F_TTT.fits, 

where V. V is the data release version, R. R and D. D are the RA 
and DEC of the tile center in degrees (J2000.0) with 1 decimal 
place, F is the filter (u, g, r, /, or ugri), and TTT is the data prod¬ 
uct type. Table [^lists the ESO product category name, file type, 

^ http://kids.strw.leidenuniv.nl/DR2 
^ http://archive.eso.org 
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value of TTT, and an example filename for each type of data 
product. 

Eor example, the KiDS-ESO-DRl r-band stacked 
image of the tile ”KIDS_48.3_-33.r' is called 
KiDS_DRl.(S)_48.3_-33. l_r_sci.fits, and the KiDS- 
ESO-DR2 multi-band source list data file corresponding to this 
tile is called KiDS_DR2.0_48.3_-33. l_ugri_src. fits. 

2.6.2. Astro-WISE archive 

All data products can also be retrieved from the As tro-WISE 
system ( [Valentijn et al.|[2007} |Begeman et al.|[MT^ , the main 
data processing and management system used for processing 
KiDS data. The pixel data and source lists are identical to the 
data stored in the ESO Science Archive Eacility, but additionally, 
the full data lineage is available in Astro-WISE. This makes this 
access route convenient for those wanting to access the various 
quality controls or further analyse, process or data-mine the full 
data set instead of particular tiles of interest. 

All data products can be accessed through the links provided 
on the KiDS website via the DBviewer interface. Downloading 
of files, viewing inspection plots, and browsing the data lineage 
is fully supported by the web interface. 

2.6.3. Synoptic table 

A third gateway to the KiDS data is the synoptic table that is in¬ 
cluded on the KiDS website. In this table quality information on 
different data products is combined, often in the form of inspec¬ 
tion figures, offering a broad overview of data quality. 


3. Observational set-up 

The KiDS survey area is split into two fields, KiDS-North and 
KiDS-South, covering a large range in right ascension so that 
observations can be made all year round. The fields, each ap¬ 
proximately 750 square degrees in size, were chosen to overlap 
with several large galaxy redshift surveys, principally SDSS, the 
2dE Galaxy Redshift Survey (2dEGRS, [Colless et al.|2001|) and 
the Galaxy And Mass Assembly (GAMA) survey ( [Driver et ak 
201 1| ). KiDS-North is completely covered by the combination 
of SDSS and the 2dEGRS, while KiDS-South corresponds to 
the 2dEGRS south Galactic cap region. Eour out of five GAMA 
fields lie within the KiDS fields. Eigure shows the outline of 
the survey fields. 

Each survey tile is observed in the u, g, r, and i bands. Ex¬ 
posure times and observing constraints for the four filters are 
designed to match the atmospheric conditions on Paranal and 
optimized for the survey’s main scientific goal of weak gravi¬ 
tational lensing. KiDS makes use of queue scheduling, allow¬ 
ing the data requiring the best conditions to be observed when¬ 
ever these conditions are met. In order to promote building up 
full wavelength coverage as quickly as possible, pointings for 
which a subset of filters has been observed receive higher prior¬ 
ity. Unfortunately, the queue scheduling system does not allow 
prioritizing of survey tiles based on the observational progress 
of neighboring tiles. As a result, the queued tiles are observed in 
a random order, resulting in the patchy on-sky distribution visi¬ 
ble in Eig.[^ The median galaxy redshift of the final survey will 
reach 0.7 and the best seeing conditions are reserved for the r- 
band, since this functions as the shape measurement band. Each 
position on the sky is visited only once in each filter, so that the 
full survey depth is reached immediately. While this precludes 
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Table 1. Data products and file types. 


Data product 

ESO product category name 

File type 


Calibrated, stacked images 

SCIENCE.IMAGE 

FITS image 

sci 

Weight maps 

ANCILLARYWEIGHTMAP 

FITS image 

wei 

Masks 

ANCILLARY.MASK 

FITS image 

msk 

Single-band source lists 

SCIENCE.SRCTBL 

Binary FITS table 

src 

Multi-band catalogue 

SCIENCE.SRCTBL 

Binary FITS table 

src 


Notes. TTT is the three character string indicating the data product type. 


variability studies, it will allow the other science projects to ben¬ 
efit from deep data from the start. During the survey the observ¬ 
ing constraints have been fine-tuned to optimize between survey 
speed and scientific return, and Table [^lists the limits employed 
for the majority of the period during which the data presented 
here were obtained. 

Since the OmegaCAM CCD mosaic consists of 32 individual 
CCDs, the sky covered by a single exposure is not contiguous 
but contains gaps. In order to fill in these gaps, KiDS tiles are 
built up from 5 dithered observations in g, r and i and 4 in u. 
The dithers form a staircase pattern with dither steps of 25" in 
X (RA) and 85" in Y (DEC), bridging the inter-CCD gaps ( [de 
I Jong et aL]|2Q13| ), see Fig. Although filled in, the gaps result 
in areas within the footprint that are covered by fewer than 5 
exposures, thus yielding slightly lower sensitivity. Due to this 
dithering strategy the final footprint of each tile is slightly larger 
than 1 square degree: 61.9 x 65.4 arcminutes in t/; 62.3 x 66.8 
arcminutes in g, r and i. Neighbouring dithered stacks have an 
overlap in RA of 5% and in DEC of 10%. The tile centers are 
based on a tiling strategy that covers the full sky efficiently for 
VST/OmegaCAM|^ The combination of the tiling and dithering 
scheme ensures that every point within the survey area is covered 
by a minimum of 3 exposures. 


4. Data processing 

The data processing pipeline used for KiDS-ESQ-DRl/2 is 


based on the Astro-WISE optical pipeline described in McFar- 
|land et aL] ( |2013[ henceforth ME 13). Below we summarize the 
processing steps and list KiDS-specific information, covering 
the KiDS process configuration and departures from the Astro- 
WISE optical pipeline. 


4.1. Image detrending 

The first processing steps are the detrending of the raw data, 
consisting of the following steps. 

Cross-talk correction. Electronic cross-talk occurs between 
several CCDs, but most strongly between CCDs #93, #94, #95 
and #96, which share the same video board (see Fig. for CCD 
numbering scheme). Cross-talk can be both positive and nega¬ 
tive, resulting in faint imprints of bright sources on neighbour¬ 
ing CCDs (Fig. [^. Although the cross-talk is generally stable, 
abrupt changes can occur during maintenance or changes to the 
instrumentation. 

A correction is made for cross-talk between CCDs #95 and 
#96, where it is strongest (up to 0.7%). Crosstalk between a pair 

^ http: //WWW .astro.rug.nl/~omegacam/dataReduction/ 
Tilingpaper.html 
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Fig. 3. Tile footprint sky coverage for observations in g, r and i. To fill 
in the gaps between the CCDs five dither steps are done, yielding a sky 
coverage pattern where the vast majority of the field-of-view is covered 
by four or five exposures. White regions are covered by five exposures 
and increasingly darker shades of gray by four, three, two and one. The 
numbers 65 to 96 indicate the CCD numbering according to the CCD 
naming scheme and the red rectangle shows how tiles are cropped while 
creating the multi-band catalogue. 


(“source” and “target”) of CCDs is determined by measuring for 
each pixel, with a value greater than 5000 ADU in the source 
CCD, the offset of the same pixel in the target CCD from the 
median value of ah pixels in the target CCD. A straight line is 
fit to the trend of this offset as function of the pixel value of the 
source CCD. See Fig. for an example. The slope of the line 
(b) is given in Table per stable period. A separate constant is 
fit to saturated pixels in the source CCD (a in the table). To cor¬ 
rect for the crosstalk between CCDs #95 and #96 the correction 
factor is applied to each pixel in the target CCD based on the 
corresponding pixel values in the source CCD: 


/;• -f d, if Ij — /sat. ? 
h "f ^^7? ^ ^sat.9 


( 6 ) 


where // and Ij are the pixel values in CCDs i and j, /' is the 
corrected pixel value in CCD i due to cross-talk from CCD j, 
and /sat. is the saturation pixel value. 

De-biasing and overscan correction. The detector bias is sub¬ 
tracted from the KiDS data in a two step procedure. First, for 
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Table 2. KiDS observing strategy: observing condition constraints and exposure times. 


Filter 

Max. lunar 
illumination 

Min. moon 
distance (deg) 

Max. seeing 
(arcsec) 

Max. airmass 

Sky transp. 

Dithers 

Total Exp. 
time (s) 

u 

0.4 

90 

1.1 

1.2 

CLEAR 

4 

1000 

g 

0.4 

80 

0.9 

1.6 

CLEAR 

5 

900 

r 

0.4 

60 

0.8 

1.3 

CLEAR 

5 

1800 

i 

1.0 

60 

1.1 

2.0 

CLEAR 

5 

1200 



Fig. 4. Example of crosstalk between CCDs 96, 95, 94 and 93 (left to right) of OmegaCAM. Both positive and negative crosstalk due to the bright 
star in CCD 94 is visible (1). Negative crosstalk on CCD 96 due to bright (saturated and unsaturated) stars on CCD 95 is indicated (2), as well as 
positive crosstalk on CCD 95 due to bright stars on CCD 96 (3). The bright ring on CCD 94 (4) is an optical "ghost" reflection. 


Table 3. Applied cross-talk coefficients. 


Period 

95 to 96 

a &(xl0-3) 

96 to 95 

a &(xl0-3) 

2011-08-01 -2011-09-17 

-210.1 

-2.504 

59.44 

0.274 

2011-09-17-2011-12-23 

-413.1 

-6.879 

234.8 

2.728 

2011-12-23 -2012-01-05 

-268.0 

-5.153 

154.3 

1.225 

2012-01-05 - 2012-07-14 

-499.9 

-7.836 

248.9 

3.110 

2012-07-14-2012-11-24 

-450.9 

-6.932 

220.7 

2.534 

2012-11-24-2013-01-09 

-493.1 

-7.231 

230.3 

2.722 

2013-01-09-2013-01-31 

-554.2 

-7.520 

211.9 

2.609 

2013-01-31 -2013-05-10 

-483.7 

-7.074 

224.7 

2.628 

2013-05-10-2013-06-24 

-479.1 

-6.979 

221.1 

2.638 

2013-06-24 - 2013-07-14 

-570.0 

-7.711 

228.9 

2.839 

2013-07-14 - 2014-01-01 

-535.6 

-7.498 

218.9 

2.701 


each science and calibration exposure the overscan is subtracted 
per row (no binning of rows). For consistency, all science and 
calibration data are reduced with this same overscan correction 
method. Second, a daily overscan-subtracted master bias, con¬ 
structed from ten bias frames and applying 3cr rejection, is sub¬ 
tracted. 


Flat-fielding. A single masterflat (per CCD and filter) was used 
for all data in the release. This is by virtue that the intrinsic pixel 
sensitivities can be considered constant to ~ 0.2% or better for g, 
r and i ( [Verdoes Kleijn et al.|2013 ). This is illustrated in Fig.|^ 
where a series of dome flat ratios is shown. The dome fiatfields 
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offer the optimal controlled experiment for assessing the pixel 
sensitivities, as the conditions under which they are observed are 
closely monitored and calibrated. More specifically, the calibra¬ 
tion unit has a special power supply which allows the ramping of 
the current when switching on and off and the stabilization of the 
current during the exposure, delivering an exposure level varia¬ 
tion less than 0.6 % over a period of 1 month ([Verdoes Kleijn 


et al.||2Ql3) ). The series of ratio plots in Fig. spans the full 


time period during the KiDS-ESO-DRl/2 observations, demon¬ 
strating that the peak-to-valley pixel-to-pixel variations vary less 
than 0.5% at any time and pixel position. 

For g, r and i the master flat is a combination of a master dome 
(for high spatial frequencies) and master twilight (=sky) flat- 
field (for low spatial frequencies). Both contributing flats are an 
average of 5 raw flat-field exposures with 3cr rejection. In u band 
only the twilight flats are used. 


Illumination correction. An illumination correction (also 
known as “photometric superflat”) is required to correct for il¬ 
lumination variations due to stray light in the flat field images, 
because of which the flat fields are not a correct representation 
of the pixel sensitivities and vignetting effects. The illumination 
correction is applied in pixel space, and only on the source fluxes 
(i.e., after background subtraction). A single illumination correc¬ 
tion image is used to correct the single master flats per Alter for 
ah data. The correction is determined from observations of sev¬ 
eral Landolt Selected Area (SA) fields ( |Landolt|1992| ) observed 
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Fig. 6. Dome fiat variation over time during the observing period of KiDS-ESO-DR2. The top panels show the 10 x 10 binned dome fiat ratios 
between a dome fiat taken on August 7 2011 and, from left to right, October 1 2011, April 1 2012, March 2 2013, and November 2 2013. The 
bottom panels show three cuts through the ratio images on top, along the lines indicated in blue on the corresponding top panels. Offsets of -0.005 
and -h0.005 are applied to the bottom and top lines, respectively, to improve legibility. 



Fig. 5. Measurement of crosstalk between CCDs #95 and #96 for one 
exposure. A straight line (green dashed) is fit to the offset from the me¬ 
dian in the target CCD (here: #95) as a function of signal in the same 
pixel in the source CCD (here: #96). Only the blue data points are used 
for this fit since at low ADU levels outliers are often present and at high 
ADU levels non-linearity affects the measurements. The slope of the 
green dashed line corresponds to the correction coefficient b. The ma¬ 
genta dotted line shows the constant a value used for saturated pixels. 


at 33 dither positions, such that the same stars are observed with 
all CCDs. After computing zero-points per CCD the residuals 
between the measured stellar magnitudes and their reference val¬ 
ues sample the illumination variations. Per CCD the illumination 
cor rection is characterized to b etter than 1%. For further details 
see |Verdoes Kleijn et al.| ( |201^ . 

De-fringing. De-fringing is only needed for the KiDS /-band. 
Analysis of nightly fringe frames showed that the pattern is 
constant in time. Therefore, a single fringe image was used 
for all KiDS-ESO-DRl images observed after 2012-01-11. For 


each science exposure this fringe image is scaled and then sub¬ 
tracted to minimize residual fringes. The general procedure, de¬ 
scribed in MF13, was modified for the KiDS data in order to 
take large-scale background fluctuations into account. Some /- 
band data contain significant background fluctuations even af¬ 
ter flatfielding due to scattered light (see Sect. |5.5| ), which can 
cause problems with the scaling of the fringe frame. Therefore, 
background-subtracted science frames were used to determine 
the scale factor. 

Pixel masking. Cosmic-rays, hot and cold pixels, and saturated 
pixels are automatically masked as described in MF13 during 
de-trending. These are included in the weight image. Additional 
automatic and ma nual masking is then applied on the coadds as 
described in Sect. 14.41 

Satellite track removal. Satellite tracks are detected by an au¬ 
tomated procedur e, working on single CCDs, that applies the 
Hough transform ( |Hough||1962[ ) to a difference image between 
exposures with the smallest dither offset. Bright stars and bright 
reflection halos are masked to limit false detections. The pix¬ 
els affected by satellite tracks are masked and included in the 
weight image. Due to the single-CCD approach small sections of 
satellite tracks in CCD corners may be missed by the algorithm. 
These remnants are included in the manual masks discussed in 
Sect.Sai 

Background subtraction. Many observations show darkened 
regions at the horizontal CCD edges where the bond wire baf¬ 
fling is placed (Iwert et al. 2006). For KiDS-ESO-DR1/2 this 
CCD-edge vignetting is corrected by performing a line-by-line 
background subtraction, implemented as a new background sub¬ 
traction method during regridding (see Sect. 3.5 of ME 13). 

The line-by-line background subtraction can be inaccurate if 
there are not enough ‘background’ pixels on a line, for example 
nearby bright stars. These regions are identified in the masks. 
Non-linearity. The linearity of the response of the OmegaCAM 
detector-amplifier chain is regularly monitored as part of the 
VST calibration procedures. No significant non-linearities are 
present (at the level of <1% over the full dynamic range of 
the system), and the pipeline currently does not include a non¬ 
linearity correction step. 
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4.2. Photometric calibration 

The steps taken to calibrate the photometry are as follows. The 
calibration described here is performed per tile and per filter and 
applied to the pixel data, which is rescaled to the fiux scale in 

Eq.[T] 


- Photometric calibration of the KiDS-ESO-DRl/2 data starts 
with individual zero-points per CCD based on SA field ob¬ 
servations (see MF13 for details of the zero-point deriva¬ 
tion). The calibration deploys a fixed aperture of 30 pix¬ 
els (~6.4" diameter) not corrected for fiux losses, and uses 
SDSS DR8 PSF magnitudes of stars in the SA fields as ref¬ 
erence. Zero-points are determined every night, except when 
no SA field observations are available, in which case default 
values, that were determined in the same way, are used. The 
total number of tiles for which these default values were used 
is 16, 13 and 14 in g, r and /, respectively, or -10% of all 
tiles. For u the nightly zero-point determinations often show 
large uncertainties and default zero-points are used more fre¬ 
quently, namely in 112 tiles. Magnitudes are expressed in AB 
in the instrumental photometric system. No colour correc¬ 
tions between the OmegaCAM and SDSS photometric sys¬ 
tem are applied. 

- Next, the photometry in the g, r and i filters is homogenized 
across CCDs and dithers for each tile independently per fil¬ 
ter. For w-band this homogenization is not applied because 
the relatively small source density often provides insufficient 
information to tie adjacent CCDs together. This global pho¬ 
tometric solution is derived and applied in three steps: 

1. From the overlapping sources across dithered exposures, 
photometric differences between the dithers (e.g. due to 
varying atmospheric extinction) are derived. 

2. Although each CCD was calibrated with its own zero- 
point, the remaining photometric differences between 
CCDs are calculated using all CCD overlaps between 
the dithered exposures. The number of sources in these 
overlaps can range between a few to several hundreds, 
depending on the filter and the size of the overlap, and 
a weighting scheme is used based on the number of 
sources. Both in steps 1 and 2, photometric offsets are ob¬ 
tained by minimizing the difference in source fluxes be- 
tween exposures and CCDs using the algorithm of |Mad-| 
[dox et al.| ( p~990| ). 

3. The offsets are applied to all CCDs with respect to the 
zero-points valid for the night, derived from the nightly 
SA field observations. 


4.3. Astrometric calibration 


A global (multi-CCD and multi-dither) astrometric c alibration is 
calculated per filter per tile. SCAMP ( |Bertin|[2006] ) is used for 
this purpose, with a polynomial degree of 2 over the whole mo¬ 
saic. The (unfiltered) 2MASS-PSC ( [Skrutskie et al.|2006| ) is used 
as astrometric reference catalogue, thus the astrometric reference 
frame used is International Celestial Reference System (ICRS). 
Using unwindowed positions the external (i.e. with respect to the 
2MASS-PSC) and internal accuracies are typically describ ed by 
a 2-D RMS of 0.3"and 0.03", respectively (see also Sect. |5.3| ). 
A more detailed description of the astrometric pipeline can be 
found in MF13. 
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4.3.1. Regridding and coadding 

SWARP ( [Bertin et al.|2002| ) is used to resample all exposures in a 
tile to the same pixel grid, with a pixel size of 0.2". After an addi¬ 
tional background subtraction step (using 3x3 filtered 128x128 
pixel blocks) the exposures are scaled to an effective zero-point 
of 0 (see also Sect. 2.4. 1| ) and coadded using a weighted mean 
stacking procedure. The applied projection method is Tangential, 
Conic-Equal-Area. Due to the individual photometric offsets ap¬ 
plied to the CCDs the gain varies slightly over the coadd. An 
average effective gain is calculated for each coadd as: 

<G>=AexpX^, (7) 

where G and S are the median gain and scale factor of the re- 
gridded CCD images and Aexp is the number of exposures. These 
average gain values are provided in the online Data tabl^ 

4.4. Masking of bright stars and defects 

Nearby saturated stars and other image defects are the main 
source of contamination in the measurement of objects. In KiDS 
coadds, these features are masked by a stand-alone program, 
PuLECENELLA (vl.O, Huang et al. in prep). Pulecenella is a novel 
procedure for automated mask creation completely independent 
of external star catalogues. An example of a mask image is 
shown in Fig. [ 7 ] 

Pulecenella detects and classifies the following types of im¬ 
age artifacts resulting from the saturated stars: 

- saturated pixels in the core and vicinity of stars, 

- spikes caused by diffraction of the mirror supports, 

- spikes caused by readout of saturated pixels, 

- “ghost” halos produced by the reflections off of optics (up to 
three wider reflection halos with spatially dependent offsets; 
these are caused by reflections of different optical elements 
in the light-path and depend on the brightness of the star). 

These features have regular shapes, and scale with the bright¬ 
ness and position of bright stars in a stable way in all images in a 
given observation band. Hence, Pulecenella is first configured 
to model the mask shapes (including the radius of saturation 
cores, the orientation of diffraction spikes, the size and offset 
of reflection halos) from some sampled saturated stars; the so 
configured analytical models are then applied to the batch mask¬ 
ing of coadds from the same band. The detection, location and 
magnitude of saturated st ars are derived from a first SExtractor 
( [Bertin & Amouts|1996| ) run over the image to be masked, with 
a band-specific configuration aiming only for the detection and 
measurement of the nearly saturated stars; the saturation pixel 
level is derived from the FITS image header. Thus, Pulecenella 
produces star masks specifically for each image, without any de¬ 
pendence on an external star catalogue. This also avoids the am¬ 
biguity in the determination of the magnitude cutoff of saturated 
stars due to the difference between the observed and external 
catalogue filters. 

The position of reflection halos is offset from the center of 
the host saturated star; the offset can be towards or outward from 
the image center, depending on the reflection components in the 
optics. For the primary reflection halo, the offset is first linearly 
modeled from several primary halos of the brightest stars in the 
image, and whether the halo mask is applied depends on the 
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Fig. 7. Examples of a Pulecenella vLO mask. Top left, thumbnail of a KiDS stacked image with several bright stars; critical areas detected by 
PuLECENELLA are overplotted with green outlines (region format). Top right: FLAG image corresponding to the image shown on the top left, with 
different colours corresponding to the pixel value with the binary flag values summed. Bottom panels: as top panels, but for a small area in a 
(different) KiDS stacked image. 


brightness of the saturated stars; this brightness level is deter¬ 
mined as the number of pixels for which the count level ex¬ 
ceeds 80% of the saturation limit and is therefore independent 
of the photometric calibration. The same method is applied for 
secondary and tertiary halos independently. The details of the 
masking method are presented in Huang et al. (2015, in prep). 


Besides the analytical star masks, Pulecenella is also config¬ 
ured to mask pixels in the empty boundary margins, CCD gaps 
or dead pixels with zero weight. These pixels are flagged as ’bad 
pixels’. Finally, bad regions which are missed by Pulecenella 
(e.g. large-scale background artifacts) but are detected by a vi¬ 
sual inspection, are manually masked using DS9 and added both 
to the region file and to the final flag image. 


As output, Pulecenella generates both an ASCII mask re¬ 
gion file which is compatible with the DS^tool and a FITS flag 
image that can be used in SExtractor. Different types of masked 
artifacts/regions are coded with the different binary values listed 
in Table 1^ In the flag image these binary values are summed. 

The flag image is used during source extraction for the 
single-band source list (see below) to flag sources whose 
isophotes overlap with the critical areas. The resulting flags are 
stored in the following two SExtractor parameters: 

- IMAELAGSJSO: sum of all mask flags encountered in the 
isophote profile, 

- NIMAELAG_ISO: number of flagged pixels entering 
IMAELAGSJSO. 


^ http://ds9.si.edu 
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Fig. 8. High-confidence star candidates and star/galaxy separation. Left: the high-confidence star candidates (red dots) are used to locate the stellar 
locus and calculate the average FWHM of the image. Right: example of star/galaxy separation; at SNR >50, the high-confidence star candidates 
(red dots) are used to define the blue line; at lower SNR, all sources with CLASS_STAR>0.8 are used; sources above the blue line are classified 
as stars. 


Table 4. Critical areas in the masks and their flag values. 


Type of area 

Flag 

Type of area 

Flag 

Readout spike 

1 

Secondary halo 

16 

Saturation core 

2 

Tertiary halo 

32 

Diffraction spike 

4 

Bad pixel 

64 

Primary halo 

8 

Manually masked 

128 


Table 5. Percentage of masked area. 


Filter 

Not masked 

Automatically 

masked 

Manually 

masked 

u 

97 

1 

2 

8 

86 

7 

7 

r 

78 

12 

10 

i 

77 

9 

14 


Table summarizes the percentages of the total area in 
KiDS-ESO-DRl/2 that are not masked, masked automatically 
by PuLECENELLA, and manually masked. Due to the lower sensi¬ 
tivity in u the number of saturated stars is much smaller, leading 
to a significantly smaller percentage of masked pixels. The high 
fraction of area that is manually masked in /-band is due to the 
high er frequency and severity of scattered light issues (see Sect. 
|5.5| ) caused by moonlight and higher sky brightness. Taking into 
account the nominal area covered by KiDS-ESO-DRl/2 the total 
unmasked area is approximately 120 square degrees. 

4.5. Source extraction and star/galaxy separation 

Single-band source lists were included in KiDS-ESO-DRl, 
while KiDS-ESO-DR2 contains both single-band source lists, as 
well as a multi-band source catalogue encompassing the com¬ 
bined area. The single-band source lists are provided per sur¬ 
vey tile, while the catalogue, although split into files correspond¬ 
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ing to single tiles, is constructed as a single catalogue. Both the 
source lists and the catalogue are intended as “general purpose” 
catalogues. 

4.5.1. Single-band source lists 

Source list extraction and star/galaxy (hereafter S/G) separation 
is achieved with an automated stand-alone procedure optimized 
for KiDS data: KiDS-CAT. This procedure, the backbone of 
which is formed by SExtractor, performs the following steps 
separately for each filter. 

1. SExtractor is run on the stacked image to measure the 
EWHM of all sources. High-confidence star candidates are 
then identified based on a number of criteria including 
signal-to-noise ratio (SNR) and ellipticity cuts (for details 
see |La Barbera et al.|2()08] ). 

2. The average PSE EWHM is calculated by applying the bi¬ 
weight location estimator to the EWHM distribution of the 
high-confidence star candidates. 

3. A second pass of SExtractor is run with SEEING_EWHM 
set to the derived average PSE EWHM. During this second 
pass the image is background-subtracted, filtered and thresh- 
olded “on the fly”. Detected sources are then de-blended, 
cleaned, photometered, and classifled. A number of SEx¬ 
tractor input parameters are set individually for each im¬ 
age (e.g., SEEING_EWHM and GAIN), while others have 
been optimized to provide the best compromise between 
completeness and spurious detections (see Data Quality sec¬ 
tion below). The detection set-up used is summarized in Ta¬ 
ble^ a full SExtractor conflguration file is available on- 
linq^ Apart from isophotal magnitudes and Kron-like ellip¬ 
tical aperture magnitudes, a large number of aperture fluxes 
are included in the source lists. This allows users to estimate 

http: //kids. strw. leidenuniv.nl/DRl/exainple_config. 
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Table 6. Detection set-up for KiDS-ESO-DRl and KiDS-ESO-DR2 single-band source lists. 


Parameter 

Value 

Description 

DETECT_THRESH 

1.5 

<sigmas> or <threshold>,<ZP> in mag/arcsec^ 

DETECT_MINAREA 

3 

minimum number of pixels above threshold 

ANALYSIS_THRESH 

1.5 

<sigmas> or <threshold>,<ZP> in mag/arcsec^ 

DEBLEND_NTHRESH 

32 

Number of deblending sub-thresholds 

DEBLEND_MINCONT 

0.001 

Minimum contrast parameter for deblending 

FILTER 

Y 

Apply Alter for detection (Y or N) 

FILTER_NAME 

default, conv 

Name of the file containing the Alter 

CLEAN 

Y 

Clean spurious detections? (Y or N) 

CLEAN_PARAM 

1.0 

Cleaning efficiency 

BACK_SIZE 

256 

Background mesh: <size> or <width>,<height> 

BACK_FILTERSIZE 

3 

Background filter: <size> or <width>,<height> 

BACKPH0T0_TYPE 

LOCAL 

can be GLOBAL or LOCAL 

BACKPH0T0 THICK 

24 

thickness of the background LOCAL annulus 


aperture corrections and total source magnitudes. All param¬ 
eters provided in the source lists are listed in the Data Format 
section below. 

4. S/G separation is performed based on the CLASS_STAR 
(star classification) and SNR (signal-to-noise ratio) param¬ 
eters provided by SExtractor and consists of the following 
steps: 

- In the SNR range where the high-confidence star candi¬ 
dates are located (the red dots in Fig. the bi-weight 
estimator is used to define their CLASS_STAR location, 
6, and its width, cr; a lower envelope of6-4cr is defined. 

- At SNR below that of the high-confidence star candi¬ 
dates, a running median CLASS_STAR value is com¬ 
puted based on all sources with CLASS_STAR >0.8. 
This running median is shifted downwards to match the 
6 - 4cr locus. The resulting curve (blue curve in Fig. 
defines the separation of stars and galaxies. 

The source magnitudes and fluxes in the final source lists 
are for zero airmass, but not corrected for Galactic foreground 
or intergalactic extinction. The result of the S/G classification is 
available in the source lists via the 2DPHOT flag. Flag values 
are: 1 (high-confldence star candidates), 2 (objects with FWHM 
smaller than stars in the stellar locus, e.g., some cosmic-rays 
and/or other unreliable sources), 4 (stars according to S/G sep¬ 
aration), and 0 otherwise (galaxies); flag values are summed, so 
2DPHOT = 5 signifles a high-confldence star candidate that is 
also above the S/G separation line. Table |A.1| lists all columns 
present in these source lists. 


employed for the source detection for the single-band source list 
(Table [^. Masking information is provided for all Alters in the 
same fashion as in the single-band source lists. Compared to the 
single-band so urce lists the number of measured parameters is 
reduced. Table |B.l[lists the columns present in these source lists. 


To account for seeing differences between Alters aperture 
corrected fluxes are provided in the catalogue. The aperture cor¬ 
rections were calculated for each Alter by comparing the aper¬ 
ture fluxes with the flux in a 30 pixel aperture, the aperture used 
for photometric calibration, and the aperture-corrected fluxes are 
included in the catalogue as separate columns. Source magni¬ 
tudes and fluxes are not corrected for Galactic foreground or in¬ 
tergalactic extinction. 


In order to prevent sources in tile overlaps to appear as mul¬ 
tiple entries in the catalogue, the survey tiles have been cropped 
and connect seamlessly to one another. This results in slightly 
shallower data along the edges of the tiles (Fig.|^, similar to the 
areas partially covered by CCD gaps. Overall, all included areas 
are covered by at least three exposures. In future releases this 
will be improved by combining information from multiple tiles 
in overlap regions. 


5. Data quality 

5.1. Intrinsic data quality 


4.5.2. Multi-band catalogue 

The multi-band catalogue delivered as part of KiDS-ESO-DR2 
is intended as a "general purpose" catalogue and relies on the 
double-image mode of SExtractor and also incorporates infor¬ 
mation obtained using the KiDS-CAT software described above. 
SExtractor is run four times for each tile, using the u, g, r and 
i KiDS-ESO-DRl/2 coadds as measurement images, to extract 
source fluxes in each of the Alters. The r-band coadd is used as 
detection image in ah runs, since it provides the highest image- 
quality in almost ah cases. For the same reason, several shape 
measurements are based only on the r-band data. In the future, 
it is foreseen that the catalogue S/G separation will make use 
of colour information and/or PSF modeling. Currently the S/G 
separation information included in the catalogue is the same as 
in the r-band source list. The detection set-up is identical to that 


In Fig. [^ the obtained seeing (FWHM), PSF ehipticity, and 
limiting magnitude distributions per Alter are shown, to illus¬ 
trate the raw data quality. The PSF ehipticity is defined here 
as 1 - BI A, where A and B are the semi-major and semi-minor 
axis, respectively. Limiting magnitudes are 5cr AB in a 2" aper¬ 
ture and determined by a At to the median SNR, estimated by 
l/MAGERR_AUTO, as function of magnitude. In case of the 
Alters observed in dark time (u, g, r) the FWHM distributions 
reflect the different observing constraints, with r-band taking the 
best conditions. Since /-band is the only Alter in which the data 
are obtained in bright time, it is observed under a large range 
of seeing conditions. Average PSF ehipticities are always small: 
<0.1 (the average is over the absolute value of ehipticity, re¬ 
gardless of the direction of ehipticity). The wide range of limit¬ 
ing magnitudes in /-band is caused by the large range in moon 
phase and thus sky brightness. 
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Fig. 9. Data quality for KiDS-ESO-DRl and KiDS-ESO-DR2. Left: average PSE size (EWHM) distributions; center: average PSF ellipticity 
distributions; right: limiting magnitude distributions (5cr AB in 2" aperture). The distributions are per filter: from top to bottom u, g, r, and 
/, respectively. The lighter portions of the histograms correspond to the 50 tiles in KiDS-ESO-DRl and the darker portions to the 98 tiles in 
KiDS-ESO-DR2. 


5.1.1. Point-spread-function 


Designed with optimal image quality over the full square degree 
field-of-view in mind, the VST/OmegaCAM system is capable 
of delivering images with an extremely uniform PSF. A typical 
example of the PSF ellipticity and size pattern for an observation 
taken with a nominal system set-up is shown in the left panel of 
Fig. [T^ Of course, as must be expected from a newly commis¬ 
sioned instrument, the set-up of the optical system is not always 
perfect and different imperfections can lead to a variety of devia¬ 
tions from a stable and round PSF. Most commonly encountered 
patterns are related to imperfect focus (increased PSF size either 
on the outside or on the inside of the FOV) and mis-alignment 
of the secondary mirror (increased PSF size and ellipticity along 
one edge of the FOV). 

To quantify the stability of the PSF we calculate a PSF size 
and ellipticity for every observation based on the PSF sizes in 32 
regions in the coadded image that correspond roughly to the 32 
CCDs. Systematic variations in ellipticity are usually less clear 
than in PSF size (see Fig. p^ , which is why the latter is used 
for monitoring PSF stability. The average PSF size in the 4 re¬ 
gions with the smallest PSF is subtracted from the average PSF 
size in the 4 regions with the largest PSF, and the result divided 
by the average PSF size of the whole image. The distribution 
per filter of this PSF size nonuniformity is shown versus the me¬ 
dian FWHM in Fig. m As the median PSF size (i.e. the see¬ 
ing) increases, the nonuniformity drops, as any differences due 
to optical imperfections are smoothed out. However, even dur¬ 
ing very good seeing conditions the PSF size variation rarely 
exceeds 25% and in most cases is around 10%. 
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The median PSF sizes for all tiles in the KiDS-North field are 
indicated in Fig.[^ Both in KiDS-North and KiDS-South there 
are no systematic gradients in PSF size over the survey area. 

5.2. Photometric quality 

The matter of photometric quality can be divided in two parts, 
namely the uniformity of the photometry within each tile, and 
the quality of the photometric calibration per tile/filter. Since the 
distribution of tiles included in KiDS-ESO-DRl/2 is not contigu¬ 
ous, with many isolated tiles, a complete photometric homoge¬ 
nization of the entire data set is impossible, and the photometric 
calibration is currently done per tile and per filter. The quality 
of the photometry will improve greatly in future releases when, 
using significant contiguous areas, a global calibration for the 
entire survey will be performed. 


5.2.1. Comparison to SDSS 


Both the internal photometric homogeneity within a coadd and 
the quality of the absolute photometric scale is assessed by com¬ 
paring the KiDS photometry to SDSS DR8 JAihara et al.|2011|), 
which is photometrically stable to ~1% (|Padmanabhan et al. 


2008| ). For this purpose, the aperture-corrected magnitudes in the 
multi-band catalogue were compared to PSF magnitudes of stars 
in SDSS DR8. Only unmasked stars with photometric uncertain¬ 
ties in both KiDS and SDSS smaller than 0.02 mag in g, r, and 
i or 0.03 mag in u were used. This comparison is only possi¬ 
ble for all tiles in the KiDS-North field, but since KiDS-South 
was calibrated in the same way as KiDS-North, we expect the 
conclusions to hold for all data. 
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KIDS 219.0 -0.5 r KIDS 131.0 1.5 r KIDS 129.4 2.5 i 



X (pixels) Y (pixels) X (pixels) Y (pixels) X (pixels) Y (pixels) 


Fig. 10. Examples of significant PSF size and ellipticity variation across three KiDS image stacks: KIDS_219.0_-0.5_r (left), KIDS_131.0_L5_r 
(center), and KIDS_129.4_2.5_i (right). In each subplot the upper panels show the FWHM in arcseconds and the lower panels the ellipticity, 
both plotted versus pixel coordinates X (left panels) and Y (right panels). Points correspond to the 500 brightest, unsaturated and unfiagged point 
sources in each field. 
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Fig. 11. PSF size nonuniformity vs. median PSF size for KiDS-ESO-DRl/2 tiles, per filter {top left: u; top right: g; bottom left: r; bottom right: i). 
The nonuniformity is calculated as the PSF size (FWHM) difference between the 4 CCDs with the biggest PSF and the 4 CCDs with the smallest 
PSF within a coadd, divided by the average PSF size. The colour scale indicates the number of tiles in each bin. 


The consistency of the photometric calibration is illustrated 
in Fig. 13 where histograms of the distributions of photometric 
offsets between KiDS and SDSS are shown for the overlapping 
tiles in KiDS-North. A systematic offset of ~ 0.02 mag is present 
in all filters, possibly due to the fact that nightly zero-points are 
determined using a fixed aperture on stars in the S A field without 
aperture correction. The scatter and occasional outliers are due 
to non-photometric conditions (during either KiDS or SA field 
observations) and, particularly in case of the w-band, use of de¬ 
fault zero-points. In Fig.f^the photometric offsets are plotted as 


function of RA and DEC, demonstrating that there are no large- 
scale gradients present. All photometric offsets determined from 


this comparison with SDSS are available in the Source catalogue 
table on the KiDS websit^El 

Figurepl]shows the residuals between KiDS and SDSS DR8 
magnitudes for one tile (KIDS_129.0_-0.5), which is a repre¬ 
sentative example. Generally speaking, the photometry in a filter 
within one tile is uniform to a few percent. The right column of 
the left panel in Fig.[^shows the distribution of the Median Ab¬ 
solute Deviation (MAD) of the stellar photometry between KiDS 
and SDSS for the tiles in KiDS-North, demonstrating the photo¬ 
metric stability within survey tiles. The relatively poor photom¬ 
etry in t/-band is due to the lack of photometric homogenization 
within a tile. 


http://kids.strw.leidenuniv.nl/DR2 
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Fig. 12. PSF sizes and photometric calibration for all tiles in KiDS-North as function of right ascension and declination. The average PSF size 
in each coadded image is indicated by the size of the circle, and the magnitude offset with respect to SDSS by its color. From top to bottom the 
panels correspond to the w, g, r, and i filters. 



Fig. 13. Comparison of KiDS aperture-corrected photometry to SDSS 
DR8 PSF photometry for stars in KiDS-North. Distributions of the me¬ 
dian offsets are shown in the left column and the distributions of Me¬ 
dian Absolute Deviations from these offsets in the right column. In both 
columns the subpanels correspond, from top to bottom, to w, g, r and /. 


5.2.2. Stellar locus 

A second quality assessment of the photometry is done by com¬ 
paring stellar photometry to empirica l stellar loci, in the f orm of 
the "principal colours" as defined by |Ivezic et al.| ( [2004| ) based 
on SDSS photometry. Only unmasked and unflagged stars with 
r<21 are used for this analysis. Although small colour terms ex- 
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Fig. 14. Photometric comparison between aperture-corrected KiDS 
photometry and SDSS DR8 PSF photometry for (unmasked) stars in 
tile KIDS_129.0_-0.5. Top: offsets versus right ascension. Bottom: off¬ 
sets versus declination. Subpanels correspond, from top to bottom, to w, 
g, r and /, respectively. Each dot corresponds to a star, with the aver¬ 
age indicated by a blue dashed line and with red dotted lines indicating 
-h0.05 and -0.05 magnitudes. 


ist between KiDS and SDSS the stellar loci based on SDSS are 
a powerful tool to verify the photometric stability over the cur¬ 
rently released tiles. Once more complete sky coverage allows an 
overall photometric calibration of the KiDS data, accurate stellar 
loci in the KiDS filters will be derived. 

Figure p3] shows the stellar locus for the tile KIDS_129.0_- 
0.5 in the three principal colour planes. These colour planes are 
different combinations of the ugri filters and denoted by v, w 
and X. The median p2 colours are calculated for each of p2s, 
p2w and p2x by choosing all stars within the indicated limits 
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Fig. 15. Stellar locus in “principal colour” space fo r tile KIDS_129.0_ - 
0.5. For the definition of the three colour planes see |Ivezic et al.| ( [2004] ). 
In SDSS filters the stellar locus should lie at p2s = p2w - p2x - 
0, indicated by the green dashed line. The dashed red line shows the 
median stellar locus colour and the solid red lines the range in the p\ 
colours used for this analysis. Only unfiagged stars with r < 21 are 
plotted here. 



and after clipping all stars more than 200 mmag away from the 
initial median p2 colour. In Fig. the distributions of the me¬ 
dian principal colours are shown, together with the typical width 
of the stellar locus, measured by the standard deviation. The nar¬ 
row distributions indicate that the colour of the stellar locus is 
typically stable to within 20 mmag. Also the stellar locus width 
is always stable to within 10 mmag. The large width of the stel¬ 
lar locus in p2x is caused by the fact that this part of the locus 
is made up of relatively faint stars with large photometric uncer¬ 
tainties. 

In this analysis three tiles stand out that have a median p2w 
colour of >100 mmag. Comparing these with the SDSS photom¬ 
etry as discussed above shows that one corresponds to the tile 
with the largest offset (-0.32 mag) in i and two correspond to 
tiles with large offsets (>0.1 mag) in r. 


5.2.3. Tile overlaps 


Finally, we analyze the tile-to-tile photometric offsets directly 
from the KiDS data, using the overlapping areas between the 
tiles where available. This is done per filter, and only stars with 
a magnitude brighter than 22 in the respective filter are matched. 
In u this results in a total number of sources per overlap region 
varying between 20 and 300, in g and r between 50 and 600 
and in i between 100 and 1,000. In some cases the PSF dete¬ 


riorates at the edge of a tile, as described Sect. |5.1.1[ and the 
aperture-corrected magnitudes may be affected. For this reason 
MAG_ISO magnitudes are used for this overlap analysis. The 
distributions of magnitude offsets in each filter, as determined 
from the tile overlaps, are shown in Fig. [T^ In the gri filters the 
offsets are typicaly 0.02 mag, and for u this increases to typi¬ 
cally 0.04 mag. Taking into account that the photometry in the 
tile edges is of slightly poorer quality than in the center, due to 
the fact that these data are less deep, these values are in good 
agreement with the photometric comparison to SDSS. 


F ig. 16. Stellar locus colour analysis using “principal colours” defined 
in |Ivezic et al.| ( |2004] ) and based on KiDS aperture-corrected stellar pho¬ 
tometry for all tiles. Distributions of the offsets from the expected stel¬ 
lar locus colour are plotted in the left column. In the right column the 
distributions of the widths of the stellar locus (measured by standard 
deviation) are shown. 
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Fig. 17. Tile-to-tile magnitude offsets in each filter measured from 
the overlap regions. In each available overlap region unfiagged stars 
were selected brighter than 22 mag in the filter considered, and their 
MAG_ISO magnitudes were compared, yielding a positive, mean off¬ 
set. 


— 1 - 1 - 1 - 1 - 1 - 1 - 

. r 

^— 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

, , , 

L. ■ 


Article number, page 15 of[^ 





















































A&A proofs: manuscript no. kids-drldr2-arxiv_v2 


[/) 

OQ 

O 




A 

A 

——I—h 


D-LLm^__ 


0.02 


0.04 


40 h 

DQ 30 

^ 20 |- 
D 

10 


60- 
i/i 50 
S 40 
O 30 
cn 20 
10 
0 


A, . 


40 

CQ 30- 
^ 20 
10 - 


25 

i/i 20 


OQ 


15- 


0.06 


O 
— 10 

5 
0 


i 


1IL_ 


22.0 


23.0 


24.0 


2D RMS (arcsec) 


98% completeness (mag) 



Fig. 18. Overview of the astrometric quality and completeness of KiDS- 
ESO-DRl and KiDS-ESO-DR2 (all data contained in the multi-band 
catalogue). Left: Histogram of the median relative astrometric offsets 
between the individual dithers and the corresponding coadd per sur¬ 
vey tile. Right: 98% completeness magnitude distributions for all tiles, 
based on the method of |Garilli et aLUlOO^ . In both columns the panels 
correspond to u, g, r, and i from top to bottom. 


Fig. 19. Spurious source contamination in the overlap between the 
KiDS multi-band catalogue and the CFHTLS-W2 field (this corre¬ 
sponds to the KiDS tiles KiDS_135.0_-1.5 and KiDS_136.0_-1.5). Top: 
spurious fraction vs. r-band magnitude (MAG_AUTO). Bottom: spuri¬ 
ous fraction vs. signal-to-noise in a 2" aperture. The black line corre¬ 
sponds to all sources, while the red line excludes sources in masked 
areas, and the blue line excludes sources in masked areas and sources 
with a non-zero SExtractor flag. 


5.3. Astrometric quality 

The accuracy of the absolute astrometry ("KiDS vs 2MASS") 
is uniform over a coadd, with typical 2-dimensional (2-D) RMS 
of 0.31" in g, r, and /, and 0.25" in u, in line with expectations 
based on the fact that the majority of reference stars used is rel¬ 
atively faint. The lower RMS in i/-band is most likely due to 
the fact that in this band on average brighter 2MASS sources 
are selected as reference sources. The accuracy of the relative 
astrometry ("KiDS vs KiDS"), measured by the 2-D positional 
residuals of sources between dithers, is also uniform across a 
single coadd. In Fig. the accuracy of this relative astrometry 
of all coadds is shown. The typical 2-D RMS is ~ 0.03" in all 
filters, but with a larger scatter in u due to the smaller number of 
available reference sources. 

5.4. Completeness and contamination 

Contamination of the multi-band catalogue by spurious sources 
was analyzed by means of a comparison of the overlap between 
KiDS and the CFHT Legacy Surve}{^ the main deeper survey 
overlapping with the current data releases (CFHTLS-W2, using 
their final data release TOOO'^^. For the analysis it is assumed 
that all KiDS sources not detected in CFHTLS-W2 are spuri¬ 
ous. Since some fraction of real sources might be absent in the 
CFHTLS catalogues, the spurious fractions derived should be 
considered upper limits. 

Figure [T^ shows the spurious fractions derived from this 
comparison as function of magnitude (r-band MAG_ISO) and 
signal-to-noise (in a 2" aperture). When all sources in the cata- 

http: //WWW .cfht.hawaii.edu/Science/CFHTLS/ 

http://terapix.iap.fr/cplt/T0007/doc/T0007-doc.html 
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logue are considered the fraction of spurious sources is estimated 
to be <5% down to a very low SNR of ~ 5 within a 2" aper¬ 
ture. Filtering sources based on masking information reduces 
this fraction to ~ 2%, demonstrating that caution is required 
when using faint sources in masked regions. When sources are 
filtered both on masking information as well as SExtractor de¬ 
tection flags, the spurious fractions drops even further to ~ 1%, 
yielding a very clean catalogue down to the detection limit. 

An internal estimate of the completeness for the KiDS data is 
provided per tile, based on the method of |Garilli et al.| ( p^99] ). It 
determines the magnitude at which objects start to be lost in the 
source list because they are below the brightness threshold in the 
detecti on cell. The implementation is similar to |La Barbera et al.| 
( [20I0[ ). Estimates of the completeness obtained by comparison 
to deeper CFHTLS-W2 data are consistent with these internally 
derived values. The distributions of the 98% completeness mag¬ 
nitudes for all tiles are shown in Fig.[^ Comparison with Fig.[^ 
shows that the 98% completeness limits are typically ~ I magni¬ 
tude brighter than the limiting magnitude for g, r and i and ~ 1.3 
magnitudes brighter in u. For the completeness of the multi-band 
catalogue the values for the r-band of each tile apply. 


5.5. Data foibles 

5.5.1. Scattered light and reflections 

Some of the main challenges in the analysis of early 
VST/OmegaCAM data are related to scattered light and reflec¬ 
tions. Due to the open structure of the telescope, light from 
sources outside the fleld-of-view often affects the observations. 
This expresses itself in a number of ways: 
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Fig. 20. Examples of known issues in the KiDS-ESO-DRl/2 VST/OmegaCAM data, highlighted by the red ellipses, a) Light patterns caused by 
reflections and scattered light of bright sources outside the EOV. b) Vignetting and scattering by CCD masks at the corners of the CCD array, c) 
Extended background structures caused by scattering of moonlight, d) Patterns caused by defective video board of CCD 82. Improved telescope 
baffles installed in early 2014 significantly reduce the occurrence of scattered light in future data releases. 


- reflections: in some cases strong reflected light patterns are 
seen in the focal plane; these are caused by light from bright 
point sources outside the fleld-of-view and can occur in all 
Alters. Some examples are shown in Fig. [20^. 

- vignetting by CCD masks: vignetting ano^cattering by the 
masks present at the corners of the focal plane array, and at 
the gaps between the rows of CCDs; this effect is particularly 
strong in i-band due to the bright observing conditions. The 
effect near the CCD gaps is largely corrected for, but in many 
cases the areas in the corners of the CCD array is strongly 
affected. Examples are shown in Fig. [2^. 

“ extended background artifacts: relatedlo the reflections men¬ 
tioned before, this is mostly seen in i-band and probably 
caused by moonlight. An example is shown in Fig. [2^. Most 
of these effects are not (yet) corrected for in the current data 
processing, but strongly affected regions are included in the 
image masks and affected sources are flagged in the source 
lists and catalogue. 

Improvements to the telescope baffles that were installed in early 
2014 significantly improve scattered light suppression. 


5.5.2. Individual CCD issues 

There are two issues related to individual CCDs that noticeably 
affect this data delivery: 


- CCD 82: this CCD suffered from random gain jumps and 
related artifacts until its video board was replaced on June 
2 2012. Artefacts as shown in Fig. [20| J are sometimes vis¬ 
ible in the image stacks due to this problem. Photometry in 
this CCD can be used due to the cross-calibration with neigh¬ 
bouring CCDs in the dithered exposures, but part of the CCD 
is lost. These features are included in the image masks and 
affected sources are flagged. 

- CCD 93: during a few nights in September 2011 (Early Sci¬ 
ence Time) one CCD was effectively dead due to a video 
cable problem. One observation included in this data deliv¬ 
ery does not include this CCD: the /-band observation of 
KIDS_341.2_-32.1. 
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6. First scientific appiications 

The Kilo-Degree Survey was designed for the central science 
case of mapping the large-scale matter distribution in the Uni¬ 
verse through weak gravitational leasing and photometric red- 
shifts. Several other science cases were also identified from the 
early design stages including studying the structure of galaxy 
halos, the evolution of galaxies and galaxy clusters, the stellar 
halo of the Milky Way, and searching for rare objects such as 
high-redshift QSOs and strong gravitational lenses. While for 
the main science goal of KiDS the full survey area is required, 
scientific analyses focusing on the other science cases are cur¬ 
rently already on-going. Below we demonstrate the quality and 
promise of KiDS by means of these on-going research efforts. 


6.1. Photometric redshifts and weak gravitational tensing 


An analysis of the weak gravitational leasing masses of galaxies 
and groups in the KiDS images is one of the main early scien- 
tific goals of the survey ( [Viola et al.|2Q15[[van Uitert et al.|2015] 
Sifon et al.|2Q15| ). Such an analysis requires two types of mea¬ 


surements from the KiDS data: galaxy shapes for shear measure¬ 
ments, and galaxy colours for photometric redshifts. The shapes 
are measured with a dedicated pipeline based on the CFHTLenS 
analy sis ( jHeymans et al.|[20T^ [Erben et ar]|2013[ [Miller et aL| 
[2013[ ) that combines information from individual exposures and 
avoids regridding of the pixels to maintain image fidelity as 
much as possible. The co lours and photometric redshifts are de¬ 
rived from the BPZ code ( [Benitez [2000t[Coe et al.[2006| ) applied 
to the output from a PSF-homogenized photometry pipeline — 
again an evolution from the CFHTLenS analysis pipeline, see 
[Hildebrand! et al.[ ( [2012[ ) — which runs on the calibrated stacked 
images released in KiDS-ESO-DRl/2. Further details of these 
dedicated analyses are presented in [Kuijken et al.[ ( [^15[ ). 

As an illustration of the quality of the photometric redshifts, 
in Fig. we show the angular cross-correlations w(6) of the 
positions of galaxies in different photometric redshift bins, on 
scales betw een 1 and 30 arcm inutes ([Erben et al.[2009{[Benjamin[ 
et al.|2010| ). We use Athena (jKilbinger et al.pOldy to calculate 


w(^) using the [Tandy & Szalay[ ( [1993[ ) estimator. Errors are ob¬ 
tained from jackknife resampling, with each pointing being a 
jackknife sub-sample. The figure shows a clear clustering signal 
of galaxies within the same redshift bin (panels on the diago¬ 
nal). Most off-diagonal panels show a smaller cross-correlation 
amplitude than the corresponding auto-correlations, and mostly 
this signal is seen only between neighbouring redshift bins, as 
would be caused by scatter of the photometric redshift estimates 
into neighbouring bins. The fact that no strong signal is seen fur¬ 
ther away from the diagonal shows that the level of catastrophic 
failures in the redshifts is low, and provides confidence in the 
photometric redshifts as well as the underlying photometry re¬ 
ported here. 


This photometric redshift-only cross-correlation check is 
complementary to the spectroscopi c redshift-photometric red¬ 
shift cross-correlations presented in [Kuijken et al.[ ( [2Q15[ ). The 
latter analysis has the advantage of utilising spectroscopic red¬ 
shifts, which provide a better representation of the ’absolute 
truth’; however, the spectroscopic redshifts only reach up to 
z ~ 0.5, so there is additional information provided by the pho¬ 
tometric redshift-only cross-correlations out to z ~ 1. We refer 
the interested reader to [Kuijken et al.[ ( [2015[ ) for details of the 
cross-correlation analysis between the photometric redshifts and 
the available spectroscopic redshifts. 
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6.2. Photometric redshifts from machine learning 


Apart from the photometric redshifts described in Sect. [6.1[ 
photometric redshifts are also derived from KiDS-ESO-DR1/2 
ugri photometry using the supervised machine learning model 
MLPQNA: a Multi-Layer Perceptron feed-forward neural net¬ 
work providing a general framework for representing nonlinear 
functional mappings between input and output variables. QNA 
stands for Quasi Newton Algorithm, a variable metric method 
used to solve optimization problems ( [Davidon|1991[ ) that, when 
implemented as the learning rule of a MLP, can be used to find 
the stationary (i.e. the zero gradient) point of the learning func¬ 
tion. The QNA implemented here is the L-BFGS algorithm from 
[Shanno[ ( [1970[ ). Supervised methods use an extensive set (the 
knowledge base or KB) of objects for which the output (in this 
case the redshift) is known a-priori to learn the mapping func¬ 
tion that transforms the input data (in this case the photometric 
quantities) into the desired output. Usually the KB is split into 
three different subsets: a training set for training the method, 
a validation set for validating the training in particular against 
overfitting, and a test set for evaluating the overall performance 
of the model ( [Cavuoti et al.[[2012l [Brescia et al.[[2013| ). In the 
method used here the validation is embedded into the training 
phase, by means of the standard leave-one-out k-fold cross vali¬ 
dation mechanism ( [Geisser|1975[ ). Performances are always de¬ 


rived blindly, i.e. using a test set formed by objects which have 
never been fed to the network during either training or valida¬ 
tion. The MLPQNA method has been successfully used in many 
experiments on different data sets, often composed through ac¬ 
curate cross-matching among public surveys (SDSS for galax- 


ies: [Brescia et al. [2014 

UKIDSS, SDSS, GALEX and WISE 

for quasars: 

Brescia et a 

l.|2013l CLASH-VLT data for galaxies: 

[Biviano et a 

L|20131(. 


The KB of spectroscopic redshifts was obtained by merg¬ 
ing the spectroscopic datasets from GAMA data release 2 ( [Liske 
[et al.|2014| ) and SDSS-III data release 9 ( [Ahn et al.|2012[ ), while 
for the KiDS photometry two different aperture magnitudes were 
adopted. The final KB includes the optical magnitudes {ugri) 
within 4"and 6"diameters and SDSS and GAMA heliocentric 
spectroscopic redshifts. GAMA redshifts come with the normal¬ 
ized quality flag NQ. 


A training set and a test set were created by splitting the KB 
into two parts of 60% and 40%, respectively. With these data 
sets two experiments were performed, one using only the GAMA 
high-quality (HQ, NQ > 2) spectroscopic redshifts and one us¬ 
ing a mix of GAMA and SDSS spectroscopic information. These 
experiments are illustrated in Fig. To quantify the quality 
of the results we use the standard normalized photometric red¬ 
shift error Aznorm defined as |zspec - ^photl/(^spec + 1). We find a 
Icr scatter in Aznorm of 0.027 and 0.031 and a fraction of catas¬ 
trophic outliers (|Aznorml > 0.15) of 0.25% and 0.39% for the 
high-quality GAMA and the GAMA-rSDSS spectroscopic data 
sets, respectively. The mixture of GAMA HQ -i- SDSS spectro¬ 
scopic data slightly extends the KB to higher redshifts, but since 
there is still very limited information at Zspec > 0.45, this does 
not significantly increase the performance at higher redshift. Fi¬ 
nally, the presence of objects at the minimum Zspec = 0 indicates 
a residual presence of stars within the sample. Further informa¬ 
tion about the experiments, results an d the produced catalo gue 
of photometric redshifts is reported in [Cavuoti et al.[p015[ ). A 
detailed comparison with redshifts from SED-fitting will be pre¬ 
sented in a forthcoming paper (Cavuoti et al. in preparation). 
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Fig. 21. Angular (cross-)correlation of galaxy positions in different photometric redshift bins in the KiDS-ESO-DR1/2 fields, based on photometric 
redshifts derived with the BPZ code. Note that the dotted lines at +0.01 demarcate a band within which the y-axis is linear. Only unmasked 
galaxies with extinction-corrected magnitudes brighter than r - 24 and IMAFLAGS_ISO=l are included. The errors are obtained from jackknife 
resampling, with each pointing being a jackknife sub-sample. The strongest signals are seen along the diagonal panels and drop off in amplitude 
for more widely separated photometric redshift bins, which gives confidence in the quality of the photometric redshifts as well as the underlying 
photometry. 


6.3. Galaxy structural parameters and scaling relations 

Galaxies are the building blocks of the Universe and within the 
framework of hierarchical structure formation they form bottom- 
up, with smaller objects forming first, then merging into massive 
structures (e.g. Lucia et al.|2006t [Trujillo et al.|2006| ). How¬ 
ever, high-mass galaxies seem to have formed most of their stars 
in earlier epochs and over a shorter time interval than the lower 
mass ones (“downsizing” scenario, [Fontanot et al.|200^ . Char¬ 


acterising the properties of the luminous matter and the way this 
has been assembled into dark-matter haloes is crucial to recon¬ 
cile theory and observations. Making use of the high spatial res¬ 
olution, depth and area coverage of KiDS we aim to measure, 
across different redshift slices: 


- total luminosity and the stellar mass through SED fitting (e.g. 

Le PHARE, [AHTouts et al.|1999l; _ 

- stellar population properties ( Bruzual & Charlot|2003| ); 
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^spec 


Fig. 22. The Zphot vs Zspec scatter plot of the blind test set KiDS-based 
optical photometric redshifts against the GAMA (black dots) and SDSS 
(red dots) spectroscopic redshifts, where the photometric redshifts have 
been obtained with the MLPQNA model. The dotted lines are the bor¬ 
ders delimiting the catastrophic outlier sub-regions (i.e. including ob¬ 
jects with \Az_norm\ > 0.15). 


- structural parameters and colour gradients. 


This will be one of the first samples of this depth over such a 
large area and wav elength coverage, allow ing us to extend pre¬ 
vious analyses (e.g. |La Barbera et al.|2010| ) to higher redshifts. 

From the KiDS-ESO-DRl/2 catalogues we select galaxies 
based on the S/G separation discussed in Sect. |4.5| with addi¬ 
tional selection criteria on size (to further reduce contamination 
by stars), and remove any flagged or masked objects. For KiDS- 
ESO-DRl/2 this results in a sample of ~ 6.5 million galaxies. 
Eor ~ 1 million of these we have measured photometric redshifts 
as discussed in Sect. |6.2| Eor the structural parameter deriva¬ 
tion we have considered only galaxies with high r-band SNR 
(SNR^ > 50) to reliably perform the surface brightness analy¬ 
sis ( |La Barbera et al.||20Q8[|2Q10| . This sub-sample consists of 


around 350,000 galaxies. The c ompleteness of the whole sample 
has been discussed in Sect. |5.4[ Eor the sample with photometric 
redshifts the 98% completeness magnitudes (all magnitudes used 
here are MAG_AUTO) are u=223, g=22.1, r=20A and /=19.7, 
which correspond to 90% completeness magnitudes of m=23.1, 
g=23.0, r=22.1 and /=21.2. The stringent cuts on SNR for the 
sample with structural parameters results in 98% completeness 
magnitudes of w=21.4, g=21.1, r=20.1, /=19.3 and 90% com¬ 
pleteness magnitudes of u=22.2, g=21.7, r=21.0, i=20.2. These 
completeness magnitudes can be converted to absolute magni¬ 
tudes sls M = m - DM(zphot) - ^corr, where DM(zphot) is the 
distance modulus based on the photometric redshift and ^corr is 
the K-correction, which is derived for two empirical models (el¬ 
liptical and Scd galaxy). Eor these calculations a standard cos¬ 
mology with Qm = 0-3, Oa = 0.7 and Hq = 75 kms“^ Mpc“^ 
is used, and the galactic foreground ext i nction is corrected based 
on the dust maps from |Schlegel et alT] ( |1998| ). Eigure [23] shows 
how the 90% completeness limits vary with absolute magnitude 
as a function of redshift. Comparing this to the M* from the lu¬ 
minosity function of low-z SDSS galaxies ( [Blanton et al.|2005] ) 
we find that M* is reached at z ~0.22, 0.31, 0.42, and 0.42 in our 
high-SNR galaxy sample and z -0.27, 0.44, 0.53, and 0.57 for 
the photo-z sample, in u, g, r and /, respectively, if the Elliptical 
model is adopted. Slightly larger redshifts are found for late-type 
galaxies (see Eig.l 
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Fig. 23. 90% completeness of the photo-z (lighter colours) and high- 
SNR (darker colours) galaxy samples in absolute magnitude vs red¬ 
shift. Solid and dashed lines are for K-corrections from Elliptical and 
Scd models. Horizontal lines indicate from the luminosity function 
analysis in jBlanton et al.| l [2b05| ). 


We make use of Galfit ( jPeng et al.|2010 ) and 2DPHOT (La 
Barbera et al.||2008| ) to fit PSE convolved Sersic profiles to the 


surface photometry and infer the structural parameters (Sersic 
index, n, effective radius Rq, axial ratio, q, disky/boxy coeffi¬ 
cient ( 24 , disk/bulge separation). In Eig.|^we show an example 
of the results obtained with 2DPHOT. The relationship between 
Rq or Sersic index and mass, luminosity and its evolution with 
redshift has been demonstrated to be a fundamental probe of 
gala xy evolution and the role of mass accretion in galaxy merg - 
ing ( [Trujillo et al.||2007] [Hilz et alT][2013[ Tortora et al.[[2014[ ). 


This characterization of 2D light profiles also allows one to de- 


termine the colour gradients ([La Barbera & de Carvalho[ 

[2009 

Tortora et al. [20lOal [Tortora et al.|201 lal[La Barbera et al.| 

2011 

Tortora & Napolitano|2012[[La Barbera et al. 

2012[), which can 

be compared to simulations (e.g Tortora et al. 

2011b Tortora et 


The resulting data set is introduced in [Tortora et al.[ ( [2015[ ) 


and applied to a first census of compact galaxies, a special class 
of objects, relic remnants of high-z red nuggets, which can pro¬ 
vide significant constraints on the galaxy merging history. The 
full analysis of the galaxy sample discussed above will be pre¬ 
sented in Napolitano et al. (in prep). 


6.4. Galaxy cluster detection 

The evolution of the number density of massive galaxy clusters 
is an important cosmological probe (e.g. [Allen et al.|20TT] ), par¬ 
ticularly at high redshifts (z > 0.5). It is therefore important to 
increase the number of known galaxy clusters. The number of 
clusters as a function of mass and redshift detected in a KiDS- 
like survey was evaluat ed using the mock catalogues derived by 
[Henriques et al.[pQ12[ HI2 hereafter), which are based on the 
semi-analytic galaxy models built by [Guo et al.[(2011[ ) for the 
Millennium Simulation (MS; [Springel et al.[[2005[ ). In particu- 
lar, these mock catalogues provide SDSS ugri photometry for 
galaxies in 24 light cones, 1.4 xl.4 deg^ each, as well as the 
M 200 mass of their parent halos. Galaxy cluste rs were identified 
in the mock catalogues following the recipe by [Milkeraitis et al.' 
( [2010| see their Table 1): cluster members were defined as those 
with the same friends-of-friends identification number and M 200 
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Fig. 24. Example of galaxy surface photometry fitting using 2DPHOT. 
Top: the thumbnail of an original KiDS r-band image with a galaxy in 
the center is shown on the left, and on the right the same thumbnail 
after subtracting the best-fit model. Bottom: normalized residuals after 
model subtraction, dljcri, are plotted as function of the distance from 
the galaxy center. Each subpanel corresponds to a different bin of the 
polar angle, 6, measured counterclock-wise from the positive x-axis. 
The residuals, normalized by the noise expected from the model in each 
pixel, CT/, are remarkably consistent with zero for all positions in the 
image, implying a reduced of about one. 


of their parent halo, and clusters with less than five members 
were rejected. Next, galaxies with gri magnitudes brighter than 
the KiDS limiting magnitudes (see Fig.|9]) were selected: we fur¬ 
ther considered only those clusters with at least 10 galaxies after 
the latter selection. 

Figure shows the mass and redshift distribution of the 
simulated clusters in H12 before (left) and after {right) the KiDS 
magnitude cuts were applied: it can be seen that a KiDS-like 
survey can probe galaxy clusters in the z ~ 0.4 - 0.8 redshift 
range, extending the cluster detection studies ba sed on SPSS 
data, which are incomp lete beyond z ~ 0.35 (e.g. [Rykoff et al.| 
|2014[[Rozo et al.|2014| ). For comparison, in the lower right panel 
we overplot the redshift distribution obtained adopting the SPSS 
limiting magnitudes (/ < 21 mag, gr < 22 mag), showing the 
improvement expected with KiDS for z > 0.4. At the moment, 


a cluster detection analysis is on-going ( 

Radovich et al.|[2015) 

using the methods described in|Bellagam 

3a et al. 

(|2011|), which 


are based on an optimal filter to find galaxy overdensities from 
the position, photometry and possibly also photometric redshift 
of the galaxies in the catalogues. Details and first results will be 
discussed in a separate paper. 


Fig. 25. Top: Mass distribution of clusters in the H12 catalogues in the 
redshift range 0 < z < 1, before (left) and after (right) the selection 
based on the KiDS gri limiting magnitudes. Bottom: Redshift distri¬ 
bution; in the lower right panel, the results obtained using the SDSS 
magnitude limits are displayed for comparison (red histogram). 


terstellar medium. The existing ensemble of high-redshift QSOs 
is dominated by luminous objects with exceptionally high ac¬ 
cretion rates (close to Eddington) and very large SMBHs (~ 10^ 
Mo). For several reasons it is needed to probe fainter QSOs, well 
below the "tip of the iceberg". This would allow to test SMBH 
early growth scenarios, which can predict a more common pop¬ 
ulation of faint QSOs ( |Costa et ar]|2014| ) with lower accretion 
rates. Furthermore, it would allow a better study of the sym¬ 
biosis in growth of SMBHs and their stellar hosts over cosmic 
time. These are hampered currently by potentially severe selec¬ 
tion biases when comparing AGN-selected sample s at high red- 
shift to the host selected samples at low redshift (|Willot et al. 
[20051 ILauer et al.|2007| ). 


For these reasons we are building up a homogeneous sam¬ 
ple of QSOs at 5.7 < z < 6.4 by combining KiDS and VIKING 
and using the /-band drop-out technique. The 9-band u through 
Ks photometry from the combined surveys goes up to ~2 mag 
deeper than SDSS, UKIDSS and the Panoramic S urvey Tele¬ 
scope & Rapid Response System 1 (Pan-STARRSl; |B~anados et 


|al.|2014j ). So far we have discovered nine such QSOs, where the 
first four are published in jVenemans et al.| ( |2015| ). KiDS is also a 
useful ingredient in the detection of very high redshift (z > 6.4) 
QSOs with VIKING. Adding the /-band KiDS data in the photo¬ 
metric selection removes -50% of these Z-band drop-out candi¬ 
dates ( jVenemans et al.|20T3| ). 


6.5. High-redshift Quasar searches 

High redshift quasars are direct probes of the Universe less than 
1 Gyr (z > 5.7) after the Big Bang. They provide fundamental 
constraints on the formation and growth of the first supermassive 
black holes (SMBHs), on early star formation, and on the chem¬ 
ical enrichment of the initially metal-free intergalactic and in- 


6.6. Strong gravitational lens searches 

Strong gravitational lensing provides the most accurate and di- 
rect probe of mass in galaxies, groups a nd clusters of galaxies 
( [Bolton et al.|2006[|Tortora et al.|2010b] ). The deep, subarcsec- 
ond seeing KiDS images are particularly suitable for a system¬ 
atic census of lenses based on the identification of arc-like struc- 
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tures around massive galaxies, galaxy groups and galaxy clus¬ 
ters. The angular size of the Einstein ring Oe can be expressed as 
function of the velocity dispersion cTy as ( [Schneider et al.|1992] ): 


6e - 1" X 


^ds 


Do 220kms“ 


( 8 ) 


where and are the angular diameter distance between the 
lens and source and to the source, respectively. For a typical 
FWHM~ 0.7" in r-band, we can expect to detect lensing arcs 
of gravitational structures with cTy > 180 kms“^ 

A first search was based on visual inspection. A sample of 
lens candidates was selected with the following two simple se¬ 
lection criteria, using the photometric redshifts described in Sect. 
|6.2| and the KiDS-ESO-DRl/2 r-band source lists. These cri¬ 
teria are aimed at maximising the lensing probability for the 
KiDS photometric sensitivities: 1) 0.1 < Zphot < 0-5 and 2) 
r < 20. Candidates are visually assessed by inspecting colour 
images and r-band i mage s where the galaxy model obtained with 
2DPHOT (see Sect. |6.3| ) is subtracted. In a control set of 600 
candidates in the area overlapping with SDSS, 18 potential lens 
candidates were identifi ed, half of which have high significance 
( jNapolitano et al.|2015| ). An example of a good lens candidate is 
shown in Fig. |26| Another approach consists of automatic selec¬ 
tion of lens candidates by filtering the KiDS-ESO-DRl/2 multi¬ 
band catalogues for massive early-type galaxies. This is done 
using colour-colour cuts and automated SED classifiers. In addi¬ 
tion, candidates of strongly lensed quasars are identified based 
on colour and morphology. Candidates from both approaches 
have been selected further via visual inspection and are being 
followed up for spectral confirmation. 

In the future, as the data volume grows, searches for strong 
lenses in large-area surveys such as KiDS will have to be per- 
formed with (semi-)automated techniques (|Alard|2006|[More etj 


al.|2012| [Gavazzi et al.|2Q14| [Joseph et al.|2014| ), as the numbers 
of typical host galaxy candidates can be of the order of thou¬ 
sands per square degree, rendering visual inspection prohibitive. 
However, current automatic lens-finding tools are not perfect and 
all automated techniques are still being tested on simulated data 
(e.g. [Metcalf & Petkova|20 14] [Gavazzi et al.|2014[ ). 


7. Summary and outlook 

The Kilo-Degree Survey is a 1500 square degree optical imaging 
survey in four filters (ugri) with the VET Survey Telescope. To¬ 
gether with its near-infrared sister survey VIKING, a nine-band 
ugriZYJHKs optical-infrared data set will be produced. While 
KiDS was primarily designed as a tomographic weak gravita¬ 
tional lensing survey, many secondary science cases are pursued. 

In this paper the first two data releases (KiDS-ESO-DRl and 
KiDS-ESO-DR2) are presented, comprising a total of 148 sur¬ 
vey tiles, or -160 square degrees. The data products of these two 
public data releases were produce d using a fine-tuned ver sion of 
the Astro- WISE optical pipeline ( [McFarland et al.|2013[ ), com¬ 
plemented with the automated Pulecenella masking software 
and the KiDS-CAT source extraction software. Data products 
include calibrated stacked images, weight maps, masks, source 
lists, and a multi-band source catalogue. Data can be accessed 
through the ESO Science Archive, the Astro-WISE system, and 
the KiDS website (see Sect. [2.6[ for the relevant links). 

The achieved data quality is very close to expectations, with 
a very small PSF ellipticity over the full FOV and typical limiting 
magnitudes (Scr AB in 2" aperture) of 24.3, 25.1, 24.9 and 23.8 
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a) b) 



Fig. 26. Example of strong lensing candidate identified in KiDS. a) 
colour image; b) r-band image; c) r-band image after subtraction of 
2DPHOT galaxy model; d) 2DPHOT galaxy model. A faint blue arc 
is visible above the central galaxy and is clearly detected after model 
subtraction. 


in u, g, r and /, respectively. As the designated lensing band, 
r-band receives the best dark seeing conditions, resulting in an 
median PSF size of -0.65", with PSF size variations over the 
FOV usually <15%. Within survey tiles, the photometry is stable 
to -2%, but since the photometric calibration is performed using 
nightly zero-points, photometric offsets between tiles are present 
due to extinction variations; this will be improved in future data 
releases. The multi-band source catalogue provided is based on 
detection in r-band and 98% complete to r ^ 24.0. 


Early weak-lensing applications ( [Sifon et al.|2015[ [Viola et[ 
[al.|2015[[van Uitert et al.|2015] ) of KiDS data rely on the KiDS- 
ESO-DRl and -DR2 data products presented here for photomet¬ 
ric redshifts, while galaxy shape measurements are based on a 
dedicated pipeline (see Erben et ah] 2013 [ [Miller et'ar[2013[ ). 
Other applications of KiDS-ESO-DRl and -DR2 data include 
galaxy structural parameter studies, cluster detection, strong 
gravitational lens detection, and high-redshift QSO searches. 
The latter has resulted in the detection of nine 5.7 < z < 6.4 
QSOs to date (Venemans et al. in prep). 

Future data releases, apart from adding more survey tiles to 
the currently covered area, are foreseen to include additional, 
value-added data products. During the first years of survey oper¬ 
ations the observational data rate has been hampered by several 
factors, causing slower progress than anticipated. Together with 
ESO the KiDS team has been working on various improvements, 
which has led to a significant increase in survey speed since early 
2014. At the current rate the survey is expected to be completed 
by 2019, but continuing efforts to enhance telescope and opera¬ 
tional efficiency are expected to lead to further improvements. 
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Appendix A: Single-band source list columns 

The following table lists the columns that are present in the single-band source lists provided in KiDS-ESO-DRl/2. Note that of the 
27 aperture flux columns only the ones for the smallest aperture (2 pixels, or 0.4" diameter) and the largest aperture (200 pixels, or 
40" diameter) are listed. Note: the label for the aperture of 28.5 pixels is FLUX_APER_28p5. 


Table A.l. Columns provided in the single-band source lists. 


Label 

Format 

Unit 

Description 

2DPHOT 

J 

Source 

classiflcation (see section on star/galaxy separation) 

XJMAGE 

E 

pixel 

Object position along x 

YJMAGE 

E 

pixel 

Object position along y 

NUMBER 

J 


Running object number 

CLASS STAR 

E 


SExtractor S/G classifler 

FLAGS 

J 


Extraction flags 

IMAFLAGSJSO 

J 


FLAG-image flags summed over the iso. profile 

NIMAFLAGJSO 

J 


Number of flagged pixels entering 

IMAFLAGSJSO 

J 


Number of flagged pixels entering IMAFLAGSJSO 

FLUX_RADIUS 

E 

pixel 

Fraction-of-light radii 

KRON_RADIUS 

E 

pixel 

Kron apertures in units of A or B 

FWHMJMAGE 

E 

pixel 

FWHM assuming a gaussian core 

ISOAREAJMAGE 

J 

pixel^ 

Isophotal area above Analysis threshold 

ELLIPTICITY 

E 


1 - B JMAGE/AJMAGE 

THETAJMAGE 

E 

deg 

Position angle (CCW/x) 

MAG_AUTO 

E 

mag 

Kron-like elliptical aperture magnitude 

MAGERR_AUTO 

E 

mag 

RMS error for AUTO magnitude 

ALPHA_J2000 

D 

deg 

Right ascension of barycenter (J2000) 

DELTA_J2000 

D 

deg 

Declination of barycenter (J2000) 

FLUX_APER_2 

E 

count 

Flux vector within circular aperture of 2 pixels 

FLUX_APER_200 

E 

count 

Flux vector within circular aperture of 200 pixels 

FLUXERR_APER_2 

E 

count 

RMS error vector for flux within aperture of 2 pixels 

FLUXERR_APER_200 

E 

count 

RMS error vector for flux within aperture of 200 pixels 

MAGJSO 

E 

mag 

Isophotal magnitude 

MAGERRJSO 

E 

mag 

RMS error for isophotal magnitude 

MAGJSOCOR 

E 

mag 

Corrected isophotal magnitude 

MAGERRJSOCOR 

E 

mag 

RMS error for corrected isophotal magnitude 

MAG_BEST 

E 

mag 

Best of MAG_AUTO and MAG JSOCOR 

MAGERR_BEST 

E 

mag 

RMS error for MAG_BEST 

BACKGROUND 

E 

count 

Background at centroid position 

THRESHOLD 

E 

count 

Detection threshold above background 

MU_THRESHOLD 

E 

arcsec"^ 

Detection threshold above background 

FLUX_MAX 

E 

count 

Peak flux above background 

MU_MAX 

E 

arcsec"^ 

Peak surface brightness above background 

ISOAREA_WORLD 

E 

deg^ 

Isophotal area above Analysis threshold 

XMINJMAGE 

J 

pixel 

Minimum x-coordinate among detected pixels 

YMINJMAGE 

J 

pixel 

Minimum y-coordinate among detected pixels 

XMAXJMAGE 

J 

pixel 

Maximum x-coordinate among detected pixels 

YMAXJMAGE 

J 

pixel 

Maximum y-coordinate among detected pixels 

X_WORLD 

D 

deg 

Baryleft position along world x axis 

Y_WORLD 

D 

deg 

B ary left position along world y axis 

XWINJMAGE 

E 

pixel 

Windowed position estimate along x 

YWINJMAGE 

E 

pixel 

Windowed position estimate along y 

X2_IMAGE 

D 

pixel^ 

Variance along x 

Y2_IMAGE 

D 

pixel^ 

Variance along y 

XYJMAGE 

D 

pixel^ 

Covariance between x and y 

X2_WORLD 

E 

deg^ 

Variance along X-WORLD (alpha) 

Y2_WORLD 

E 

deg^ 

Variance along Y-WORLD (delta) 

XY_WORLD 

E 

deg^ 

Covariance between X-WORLD and Y-WORLD 

CXXJMAGE 

E 

pixeU^ 

Cxx object ellipse parameter 

CYYJMAGE 

E 

pixeU^ 

Cyy object ellipse parameter 
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Table A. 1 - continued from previous page 


Label 

Eormat 

Unit 

Description 

CXYJMAGE 

E 

pixeL^ 

Cxy object ellipse parameter 

CXX_WORLD 

E 

deg“^ 

Cxx object ellipse parameter (WORLD units) 

CYY.WORLD 

E 

deg"^ 

Cyy object ellipse parameter (WORLD units) 

CXY.WORLD 

E 

deg"^ 

Cxy object ellipse parameter (WORLD units) 

AJMAGE 

D 

pixel 

Profile RMS along major axis 

BJMAGE 

D 

pixel 

Profile RMS along minor axis 

A_WORLD 

E 

deg 

Profile RMS along major axis (WORLD units) 

B_WORLD 

E 

deg 

Profile RMS along minor axis (WORLD units) 

THETA_WORLD 

E 

deg 

Position angle (CCW/world-x) 

THETA_J2000 

E 

deg 

Position angle (east of north) (J2000) 

ELONGATION 

E 

deg 

A_IMAGE/B_IMAGE 

ERRX2_IMAGE 

E 

pixel^ 

Variance of position along x 

ERRY2_IMAGE 

E 

pixel^ 

Variance of position along y 

ERRXYJMAGE 

E 

pixel^ 

Covariance of position between x and y 

ERRX2_WORLD 

E 

deg2 

Variance of position along X-WORLD (alpha) 

ERRY2_WORLD 

E 

deg2 

Variance of position along Y-WORLD (delta) 

ERRXY_WORLD 

E 

deg^ 

Covariance of position X-WORLD/Y-WORLD 

ERRCXXJMAGE 

E 

pixeL^ 

Cxx error ellipse parameter 

ERRCYYJMAGE 

E 

pixeL^ 

Cyy error ellipse parameter 

ERRCXYJMAGE 

E 

pixeL^ 

Cxy error ellipse parameter 

ERRCXX_WORLD 

E 

deg“2 

Cxx error ellipse parameter (WORLD units) 

ERRCYY_WORLD 

E 

deg“^ 

Cyy error ellipse parameter (WORLD units) 

ERRCXY_WORLD 

E 

deg“^ 

Cxy error ellipse parameter (WORLD units) 

ERRAJMAGE 

E 

pixel 

RMS position error along major axis 

ERRBJMAGE 

E 

pixel 

RMS position error along minor axis 

ERRA_WORLD 

E 

deg 

World RMS position error along major axis 

ERRB_WORLD 

E 

deg 

World RMS position error along minor axis 

ERRTHETAJMAGE 

E 

deg 

Error ellipse position angle (CCW/x) 

ERRTHETA_WORLD 

E 

deg 

Error ellipse position angle (CCW/world-x) 

ERRTHETA_J2000 

E 

deg 

J2000 error ellipse pos. angle (east of north) 

EWHM_WORLD 

E 

deg 

EWHM assuming a gaussian core 

ISOO 

J 

pixel^ 

Isophotal area at level 0 

ISOl 

J 

pixel^ 

Isophotal area at level 1 

IS02 

J 

pixel^ 

Isophotal area at level 2 

ISOS 

J 

pixel^ 

Isophotal area at level 3 

IS04 

J 

pixel^ 

Isophotal area at level 4 

1805 

J 

pixel^ 

Isophotal area at level 5 

IS06 

J 

pixel^ 

Isophotal area at level 6 

ISO? 

J 

pixel^ 

Isophotal area at level 7 

SLID 

K 


Source list ID 

SID 

K 


Source ID within the source list 

HTM 

K 


Hierarchical Triangular Mesh (level 25) 

ELAG 

K 


Not used 
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Appendix B: Multi-band catalogue 

The following table lists the columns that are present in the multi-band catalog provided in KiDS-ESO-DR2. 
Table B.l. Columns provided in the multi-band catalogue. 


Label 

Eormat 

Unit 

Description 

Measurements based on r-band detection image 

ID 

23A 


Source identifier 

RAJ2000 

D 

deg 

Right ascension (J2000) 

DECJ2000 

D 

deg 

Declination (J2000) 

SG2DPHOT 

K 


Source classification 

A 

D 

pixel 

Linear semi major axis 

B 

D 

pixel 

Linear semi minor axis 

CLASS_STAR 

E 


SExtractor star/galaxy classifier 

KRON RADIUS 

E 

pixel 

Kron-radius used for MAG_AUTO 

POSANG 

E 

deg 

Position angle 

SEQNR 

K 


r-band sequence number 

Measurements provided for each filter 

PLUXERR_APER_100_<filter> 

E 

count 

fiux error in 100 pixel aperture 

PLUXERR_APER_10_<filter> 

E 

count 

flux error in 10 pixel aperture 

PLUXERR_APER_14_<filter> 

E 

count 

fiux error in 14 pixel aperture 

PLUXERR_APER_25_<filter> 

E 

count 

fiux error in 25 pixel aperture 

PLUXERR_APER_40_<filter> 

E 

count 

fiux error in 40 pixel aperture 

PLUXERR_APER_4_<filter> 

E 

count 

fiux error in 4 pixel aperture 

PLUXERR_APER_6_<filter> 

E 

count 

fiux error in 6 pixel aperture 

PLUXERR_APERCOR_l 00_<filter> 

E 

count 

corrected fiux error in 100 pixel aperture 

PLUXERR_APERCOR_l 0_<filter> 

E 

count 

corrected fiux error in 10 pixel aperture 

PLUXERR_APERCOR_14_<filter> 

E 

count 

corrected fiux error in 14 pixel aperture 

PLUXERR_APERCOR_25_<filter> 

E 

count 

corrected fiux error in 25 pixel aperture 

PLUXERR_APERCOR_40_<filter> 

E 

count 

corrected fiux error in 40 pixel aperture 

PLUXERR_APERCOR_4_<filter> 

E 

count 

corrected fiux error in 4 pixel aperture 

PLUXERR_APERCOR_6_<filter> 

E 

count 

corrected fiux error in 6 pixel aperture 

PLUX_APER_100_<filter> 

E 

count 

fiux in 100 pixel aperture 

PLUX_APER_10_<filter> 

E 

count 

fiux in 10 pixel aperture 

PLUX_APER_14_<filter> 

E 

count 

fiux in 14 pixel aperture 

PLUX_APER_25_<filter> 

E 

count 

fiux in 25 pixel aperture 

PLUX_APER_40_<filter> 

E 

count 

fiux in 40 pixel aperture 

PLUX_APER_4_<filter> 

E 

count 

fiux in 4 pixel aperture 

PLUX_APER_6_<filter> 

E 

count 

fiux in 6 pixel aperture 

PLUX_APERCOR_100_<filter> 

E 

count 

corrected fiux in 100 pixel aperture 

PLUX_APERCOR_10_<filter> 

E 

count 

corrected fiux in 10 pixel aperture 

PLUX_APERCOR_14_<filter> 

E 

count 

corrected fiux in 14 pixel aperture 

PLUX_APERCOR_25_<filter> 

E 

count 

corrected fiux in 25 pixel aperture 

PLUX_APERCOR_40_<filter> 

E 

count 

corrected fiux in 40 pixel aperture 

PLUX_APERCOR_4_<filter> 

E 

count 

corrected fiux in 4 pixel aperture 

PLUX_APERCOR_6_<filter> 

E 

count 

corrected fiux in 6 pixel aperture 

PLUX_RADIUS_<filter> 

E 

pixel 

SExtractor ELUX_RADIUS 

PWHM_IMAGE_<filter> 

E 

pixel 

SExtractor EWHMJMAGE 

PLAG_<filter> 

J 


SExtractor extraction fiag 

IMAPLAGS_ISO_<filter> 

J 


Mask fiag 

MAGERR_AUTO_<filter> 

E 

mag 

RMS error for MAG_AUTO 

MAGERR_ISO_<filter> 

E 

mag 

RMS error for MAG_ISO 

MAG_AUTO_<filter> 

E 

mag 

Kron-like elliptical aperture magnitude 

MAG_ISO_<filter> 

E 

mag 

Isophotal magnitude 

NIMAPLAGS_ISO_<filter> 

J 


Number of masked pixels entering IMAELAGS_ISO 

ISOAREA_IMAGE_<filter> 

J 

pixel^ 

Isophotal aperture 

XPOS_<filter> 

E 

pixel 

X pixel position <filter> coadd 

YPOS<filter> 

E 

pixel 

Y pixel position <filter> coadd 
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